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Transition  of  the  1993  Version  of  the 
Aeroorediction  Code  fAP93) 


The  AP93  code  has  recently  been  completed  and  is  now 
available  for  transition  to  requesting  users.  Users  of  the  1981 
version  of  the  code  included  approximately  60  DOD  contractors,  DOD 
or  other  government  agencies,  universities  and  several  foreign 
countries.  A  small  charge  was  incurred  by  the  requester  to 
alleviate  administrative  costs  of  the  transition.  In  an  effort  to 
promote  dual  use  technology  between  government  and  industry, 
transition  of  the  AP93  will  be  done  free  of  charge  to  requesting 
users.  This  includes  either  a  9  track  tape  or  disk  along  with  a 
users  guide.  For  DOD  contractors  or  other  government  agencies,  the 
attached  form  must  be  filled  out  completely  and  signed  by 
appropriate  authorities.  For  foreign  countries  or  companies  owned 
by  foreign  countries,  the  request  must  go  from  the  foreign  embassy 
to  the  U.S.  Embassy  for  approval  since  the  code  is  for  limited 
distribution.  The  U.S.  Embassy  will  then  forward  the  request  to 
the  Office  of  Naval  Technology  for  implementation  at  NSWCDD.  Since 
we  are  sending  out  a  copy  of  the  new  code  free  of  charge,  we  will 
not  be  able  to  provide  the  level  of  consulting  that  we  would  like 
to.  However,  we  do  request  several  things  from  the  users.  First, 
if  you  find  an  error  in  the  code,  we  would  appreciate  you  bringing 
it  to  either  my  attention  at  (703)  663-8141,  to  Mr.  Tom  Hymer  at 
(703)  663-4164,  or  to  Dr.  Roy  Mclnville  at  (703)  663-4760. 
Secondly,  if  you  find  areas  of  success  or  areas  of  weakness,  we 
would  appreciate  hearing  from  you.  We  use  information  such  as  this 
to  help  guide  any  future  modifications.  Thirdly,  if  you  find  areas 
where  additional  code  capability  would  be  desirable,  again  we  would 
like  to  hear  from  you.  In  all  cases,  a  written  letter  is  preferred 
but  a  phone  call  is  better  than  no  communication  at  all.  Also,  if 
you  use  the  code  for  the  purpose  of  comparing  the  results  to  other 
computations  within  the  external  literature,  we  would  like  you  to 
use  the  terminology  AP93.  If  you  modify  the  AP93  to  suit  your  own 
needs  you  could  refer  to  the  code  as  modified  AP93.  Finally,  we 
request  that  you  not  send  any  copies  of  the  code  to  other  agencies 
or  even  other  activities  in  your  company  if  it  is  at  a  different 
location.  Please  refer  them  to  us  so  that  we  can  keep  track  of 
copies  of  the  code. 


DR.  FRANK  MOORE 
Dahlgren  Division 
Naval  Surface  Warfare  Center 
Dahlgren,  VA  22448 


AP93  Transition  Request 


NAME:  _ 

COMPANY  NAME: 
ADDRESS : 


PHONE  NUMBER: 


GOVT  CONTRACT  NUMBER: 


MEANS  OF  TRANSITION:  9  Track  Tape  _  5. 25"  Diskette 

3.5”  Diskette  E-Mail 


Please  supply  E-Mail  address  if 

E-Mail  is  desired _ 

We  agree  to  not  give  a  copy  of  the  AP93  to  any  other  agency.  We 
will  refer  them  to  NSWCDD  at  the  address  below  for  a  copy. 

Your  Signature  _ 

Supervisor  Signature  _ 

Please  send  this  form  completely  filled  out  to: 

Commander 

Attn  Dr.  Frank  Moore  (G04) 

Dahlgren  Division 

Naval  Surface  Warfare  Center 

Dahlgren,  VA  22448 

We  will  try  to  get  you  a  copy  of  the  AP93  along  with  a  users  guide 
within  a  month. 


NSWCDD/TR-93/91 


IMPROVED  AEROPREDICTION  CODE: 
PART  l-SUMMARY  OF  NEW  METHODS 
AND  COMPARISON  WITH  EXPERIMENT 


BY  FRANK  G.  MOORE  THOMAS  C.  HYMER  ROY  M.  MCINVILLE 
WEAPONS  SYSTEMS  DEPARTMENT 


MAY  1993 


Approved  for  public  release,  distribution  is  unlimited 


NAVAL  SURFACE  WARFARE  CENTER 
DAHLGREN  DIVISION 
Dahlgren,  Virginia  22448-5000 


NSWCDD/TR-93/91 


FOREWORD 


The  effort  described  in  this  report  is  the  culmination  of  a  three-and-half-year 
effort  to  upgrade  the  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD) 
aeroprediction  code  (AP93)  to  meet  the  needs  that  have  arisen  over  the  past  10  years. 
These  needs  included  high  Mach  numbers,  accurate  nonlinear  aerodynamics  for  low- 
aspect  ratio  missile  configurations,  and  improved  base-drag  prediction.  This  report 
summarizes  all  the  new  technology  developed  and  shows  example  comparisons  of  the 
new  AP93  code  to  the  former  1981  version.  A  user’s  guide  and  program  listing  will  be 
included  in  Part  II  of  this  report,  which  will  have  a  more  limited  distribution. 

Other  personnel  who  participated  in  this  effort  are  Fred  DeJarnette  of  North 
Carolina  State  University,  Frank  Baltakis  of  Advanced  Technology  Associates,  Floyd 
Wilcox  of  the  National  Aeronautics  and  Space  Administration/Langley  Research 
Center,  and  Mike  Armistead  ,  Steve  Rowles,  and  Leroy  Devan  from  NSWCDD. 
Appreciation  is  expressed  to  each  of  these  individuals  for  their  roles,  which  have  been 
documented  in  previous  technical  reports. 

The  work  described  in  this  report  was  supported  through  the  Office  of  Naval 
Research  (Dave  Siegel)  and,  more  specifically,  the  Surface-Launched  Weapons 
Technology  Block  Program  managed  at  NSWCDD  by  Robin  Staton.  The  authors 
express  their  appreciation  to  these  individuals  for  their  support  in  this  work. 


Approved  by: 

DAVID  S.  MALYEVAC,  Deputy  Department  Head 
Weapons  Systems  Department 
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ABSTRACT 

A  new  and  improved  version  of  the  Naval  Surface  Warfare  Center,  Dahlgren 
Division  aeroprediction  code  (AP93)  has  been  developed.  The  new  code  contains  new 
technology  that  allows  planar  aerodynamics  of  axisymmetric  solid  rocket-type 
weapons  to  be  computed  with  engineering  accuracy  over  the  entire  Mach  number 
range  and  for  angles  of  attack  to  30  deg.  New  technology  developed  and  included  in 
the  AP93  includes 

•  A  new  engineering  method  to  compute  aeroheating  information  at  a 
high  Mach  number 

•  Extension  of  the  second-order  shock-expansion  theory  to  include 
real-gas  effects,  including  several  new  pressure  prediction  tech¬ 
niques 

•  An  improved  body-alone  nonlinear  normal-force  method 

•  New  methods  for  computing  nonlinear  aerodynamics  of  wing  alone, 
wing  body,  and  body  wing  due  to  angle  of  attack,  and  wing  body  due 
to  control  deflection 

•  A  new  base-drag  database  and  improved  empirical  base-drag 
estimation  technique 

Comparison  of  the  AP93  code  to  the  former  aeroprediction  code  (AP81)  and 
experimental  data  on  many  configurations  and  test  conditions  showed  the  following: 

•  The  AP93  code  reduces  the  normal-force  and  center-of-pressure 
errors  of  the  AP81  code  considerably  for  most  configurations. 

•  The  AP93  code  is  more  robust  in  terms  of  accuracy  over  a  broad 
range  of  Mach  numbers,  angles  of  attack,  and  configuration 
geometries. 

•  The  AP93  code  gave  slightly  improved  axial-force  coefficients  on 
average. 

•  Computational  time  and  ease  of  use  were  about  the  same. 

Comparison  of  the  AP93  code  to  available  computations  of  other  state-of-the-art  codes 
shows  the  AP93  to  be  as  good  or  superior  to  these  codes  for  planar  aerodynamics. 

The  new  code  and  associated  technology  will  be  available  for  transition  to 
legitimate  requesting  users  by  September  1993  at  no  charge  to  the  user. 
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1.0  INTRODUCTION  AND  BACKGROUND 


For  the  past  20  years,  the  Naval  Surface  Warfare  Center,  Dahlgren  Division 
(NSWCDD)  has  been  involved  in  developing  codes  to  calculate  aerodynamics  on 
tactical  weapons.  These  codes  have  attempted  to  meet  the  changing  needs  of  the 
Tactical  Weapons  Community  and  keep  pace  with  aerodynamic  requirements.  A 
recent  effort1  was  undertaken  to  look  at  where  we  have  been,  where  we  are,  and 
where  we  need  to  go  in  the  future  with  respect  to  aerodynamic  codes.  Ore  of  the 
primary  needs  identified  iu  Reference  1  was  an  upgrade  of  the  NSWCDD  aero- 
prediction  code  (APC)  to  allow  Mach  numbers  up  to  20  (including  the  effects  of  real 
gases; ,  improved  lift  prediction  with  particular  emphasis  on  low-aspect  ratio  lifting 
surfaces,  and  improved  base-drag  prediction.  All  three  of  these  efforts  were  under¬ 
taken  and  have  now  been  completed.  In  addition,  other  minor  modifications  to  the 
1981  version  of  the  APC  have  been  made,  including  options  for  boundary  layer 
transition  and  generating  aerodynamic  heat  transfer  coefficients.  This  report  serves 
as  a  summary  theoretical  document  of  the  new  APC;  i.e.,  aeroprediction  1993  (AP93). 
The  new  technology  developed  will  be  summarized,  and  the  technology  available  in 
the  1981  version  of  the  APC;  i.e.,  aeroprediction  1981  (AP81),  that  remains  in  the 
AP93  will  also  be  discussed.  Finally,  a  new  technology,  a  new  nonlinear  wing-bodv 
interference  factor  for  control  deflections,  until  now  not  documented,  will  be 
discussed  in  more  detail.  In  addition,  come  other  minor  modifications  to  the  new 
nonlinear  theory  will  be  introduced. 

The  four  former  versions  of  APC  were  documented28  and  transitioned  to  users 
in  1972,  1974,  1977,  and  1981.  Each  of  these  versions  attempted  to  meet  the 
requirements  as  seen  by  the  tactical  weapons  community.  The  first  version2  was  for 
general-shaped  bodies  alone.  It  was  the  first  such  code  known  that  combined  a  good 
mix  of  accuracy  in  aerodynamic  computations,  ease  of  use  and  computational  time.  It 
is  believed  that  this  mix  led  to  the  code’s  initial  popularity  and  requests  for  additional 
capability.  In  1974,3,4  the  code  was  extended  to  allow  up  to  two  sets  of  lifting  surfa-es 
in  the  computational  process.  In  1977,5,6  dynamic  aerodynamic  derivatives  were 
added  to  the  code’s  capability.  Finally,  the  last  version  of  the  code7-8  extended  the 
Mach  number  range  up  to  eight  and  added  high  angle-of-attack  capability  for  a 
narrow  range  of  configurations. 

Over  the  past  10  years,  the  AP81  has  been  used  to  compile  aerodynamics  on 
configurations  at  conditions  where  the  accuracy  is  not  good.  This  includes  angles  of 
attack  greater  than  about  15  deg  on  missiles  with  two  sets  of  lifting  surfaces  and 
Mach  numbers  greater  than  eight  where  real-gas  effects  become  important  in  the 
aerothermal  environment  Furthermore,  the  base-drag  estimation,  while  including 
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angle  of  attack  and  first-order  fin  effects,  needed  additional  wind  tunnel  data  for 
more  accurate  computations.  As  a  result  of  these  known  shortcomings,  a  desire  to 
still  use  the  AFC  and  the  fact  that  there  is  no  accurate  engineering  code  available  to 
accomplish  the  objectives  of  angle  of  attack,  Mach  number,  and  base-^rag  prediction 
capability,  an  effort  was  begun  in  1990  to  extend  the  APC  to  meet  these  require¬ 
ments. 


The  extension  of  the  code  to  hypersonic  Mach  numbers,  including  real-gas 
effects,  was  completed  and  documented. 9.10  In  the  development  of  this  new  tech¬ 
nology,  three  new  pressure  prediction  methods  were  derived  along  with  a  method  for 
accurately  estimating  inviscid  surface  temperatures.  The  inviscid  surface  tempera¬ 
tures  were  then  used  to  develop  an  engineering  method  for  boundary  layer  heating. H 
The  new  method  thus  provides  engineering  estimates  of  adiabatic  wall  temperature 
and  heat  transfer  coetTicients  for  configurations  flying  at  hypersonic  Mach  numbers 
and  at  angle  of  attack. 

New  methods  for  nonlinear  wing-alone,  wing-body,  and  body-wing  normal  force 
were  recently  developed. 12.13  In  addition,  an  improved  body-alone  nonlinear  normal 
force  prediction  capability  was  developed.  Improved  center-of-pressure  estimates 
were  added  that  gave  improved  pitching  moments.  The  code’s  angle-of-attack  range 
is  now  0  to  30  deg  in  the  planar  or  U|>  =  0 )  roll  position. 

Wind  tunnel  tests  were  conducted  to  measure  base  pressure  as  a  function  of 
Mach  number,  angle  of  attack,  fin  thickness,  fin  location,  and  fin  deflection  Using 
this  data,  an  improved  empirical  base-drag  prediction  methodology  has  been 
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A  new  nonlinear  method  for  wing-body  interference  due  to  control  deflection 
has  recently  been  developed.  This  new  technology,  along  with  options  for  skin 
friction  drag  computation,  will  be  discussed  in  this  report. 

Thus,  a  summary  of  the  new  technology  developed,  a  review  of  the  technology 
used  in  the  APB  I  that  will  remain  in  AP93,  along  with  results  and  discussion  on 
several  missile  configurations  will  be  given  in  this  report.  Where  results  of  other 
state-of-the-art  aerodynamics  codes  are  available,  these  results  will  also  be  presented 
for  comparison  purposes.  The  second  volume  of  this  report,  which  has  a  limited 
distribution  statement,  contains  the  computer  code  listing,  along  with  guidelines  on 
how  to  use  the  code. 
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2.0  SUMMARY  OF  NEW  METHODS  IN  AP93 


2.1  HYPERSONIC  MACH  NUMBER  METHODOLOGY 


The  hypersonic  Mach  number  methodology,  which  in  the  AP93  is  taken  to  be 
M<»5:6,  is  based  on  the  use  of  the  second-order  shock-expansion  theory  (SOSET)  with 
inclusion  of  real-gas  effects.  The  details  of  the  new  technology  are  reported  in 
Reference  9  and  summarized  in  Reference  10.  An  even  briefer  summary  follows. 


SOSET  defines  the  pressure  at  any  point  on  a  pointed  body  of  revolution  as 


P  =  Pc  -  (Pc  ~  P^e  n  1 

where  n  is  an  exponential  decay  term  given  by 

(f)2<8-82>  <21 

Refer  to  Figure  1  for  the  nomenclature. 


FIGURE  1.  FLOW  ABOUT  A  FRUSTUM  ELEMENT 


X 
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Here  r]  is  positive.  If  conditions  are  such  that  q  becomes  negative,  then 
Equation  1  does  not  satisfy  the  boundary  conditions  of  p  =  pc  as  s  — ►  <t>  and  therefore  q 
must  be  defined  by  an  equation  different  from  Equation  2.  A  value  of  q  =  0  causes  the 
SOSET  to  revert  back  to  the  generalized  shock -expansion  theory  (GSETi;  a  value  of  q 
between  q  =  0  and  q  =  o°  allows  a  blend  of  these  two  theories.  To  extend  SOSET  to  real 
gases  requires  several  things:  (Da  cone  solution  for  real  gases  <pc);  <2)  a  Prandtl- 
Mever  Expansion  (PME)  for  real  gases  ( p*2 );  <3»  a  derivation  of  a  new  pressure 
derivative  (dp/ds)g,  where  the  perfect-gas  assumption  has  not  been  made;  and  (4)  a 
way  to  compute  temperature  given  values  of  pressure.  After  the  real-gas  pressure 
derivative  (dp/<?s)2  was  derived  and  checked,  it  was  found  that  iftp.’cisD  became 
negative  for  many  cases  causing  one  to  choose  between  the  GSET  (q-0)  and  the 
tangent  cone  theory  lq  =  <»).  In  comparisons  of  the  pressure  prediction  to  full  Euler 
computations,  it  was  found  that  a  better  way  to  implement  the  shock  expansion 
theory  for  Msfi  was  to  redefine  Equation  1  as 

p  =  Pc  -  <P.  -  p2»  rj !  (3» 

with  qi  being  an  input  parameter  chosen  by  the  user.  It  was  found  that  a  value  of 
qi  =  0  gave  slightly  better  pressure  predictions  for  slightly  blunt  configurations, 
whereas  a  value  of  qi  =  1  gave  better  accuracy  where  bluntness  was  large.  The  final 
implementation  of  SOSET  in  AP93  is  thus  Equation  3  allowing  qi  as  an  input,  pc  is 
the  real-gas  tanget  cone  pressure,  and  P2  is  the  real-gas  value  of  pressure  computed 
from  a  Prandtl-Meyer  expansion. 

To  compute  inviscid  temperatures  (and  other  properties)  along  the  surface  of  a 
pointed  or  blunt  body,  use  is  made  of  the  constancy  of  entropy  along  the  surface  for 
perfect,  frozen,  or  equilibrium  chemically  reacting  flows.  Knowing  the  value  of 
entropy  and  pressure  from  the  pointed  cone  solution  of  Reference  16  or  the  normal 
shock  solution  of  Reference  17  for  a  blunt  body,  one  can  then  use  the  thermofit 
equations  of  References  18  and  19  to  determine  other  properties;  i.e.. 


T  =  Tip, Si 

(4a) 

p  =  p  (p,S) 

(4b) 

a  =  a  ip.S) 

(4c) 

e  —  e(p.S) 

(4d) 

The  remaining  properties  at  the  body  surface  can  be  found  from 
dynamic  relationships;  i.e.. 

standard  thermo- 

h  = e  +  pp 

(5a) 

H  —  I  )t  -  constant 

O  \  Y  _  1  )  0r 

*  r 

(5b) 
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V  =  v/2(H  -  h) 

O 


<5c) 


M  =  V/a 


(5d) 


Y  = 


(5e) 


Z  = 


P 

pRT 


(50 


In  the  process  of  computing  surface  properties,  three  new  pressure  prediction 
methods  were  derived.  The  first  of  these  was  to  give  an  improved  pressure  coefficient 
prediction  on  the  blunt  nose  of  a  missile  configuration  over  that  provided  by  the 
modified  Newtonian  theory  (MNT).  If  the  pressure  coefficient  of  MNT  is  defined  as 


=  C  sin25 

PQ  eq 


(6) 


then  the  new  pressure  on  the  blunt  nose  part  of  a  missile  is  given  by 


C 

p 


P'MNT 


-  AC 

p 


(7) 


ACp  of  Equation  7  is  defined  by 


ACp  =  kcosm(8eq)lcos5eq  -  cos(6eq)m! 


(8) 


where(8eq)m  =  25.95  deg,  m  =  2.78,  and 


k  =  2.416C  +  4.606 

P„ 


0.1507  C2  +  C 

p0  M2  po 


L’2 


Figure  2  shows  the  results  of  the  improved  modified  Newtonian  theory  (IMNT) 
of  Equations  7  and  8,  compared  to  Equation  6  alone  and  a  full  numerical  solution  of 
the  Euler  equations20  for  a  hemispherical  forebody  at  M*  =  10.  The  IMNT  gives  up  to 
7-percent  improvement  in  pressure  compared  to  the  MNT.  Even  past  the  match  point 
(5  eq  <  25.95  deg),  the  IMNT  gives  good  agreement  with  the  numerical  solution  down 
to  8eq  values  of  10  deg.  This  level  of  accuracy  in  pressure  prediction  will  also 
translate  into  more  accurate  drag  computations,  particularly  on  bodies  with  large 
bluntness. 
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<5Eq  (deg.) 


FIGURE  2.  SURFACE  PRESSURE  DISTRIBUTION  OVER 
A  HEMISPHERICAL  FOREBODY  AT  Mx=  10 


The  other  two  pressure  prediction  formulas  have  to  do  with  calculating  the 
pressure  on  a  point  behind  the  blunt  nose  portion  of  the  body  but  at  an  angle  of 
attack.  These  are 

C  <.ocT<j>)  =T  0  -  (2«)sin(29)cos(4>)  +  {Feos20)<*2  +  (4/3  sin  (20)  cos  ((JO)® 3  ^ 

P  P<r=0 


where 


F  =  (2-  -  |(1 


sin2<{> 


and 


Cp  («,4>) 


(2«)sin(20)cos(4») 

3 


(10) 


Equation  9  is  used  for  pointed  body  configurations,  as  well  as  for  blunt  body 
configurations  in  the  windward  plane  area  (60°<$^  180°).  Equation  10  is  used  in  the 
leeward  plane  (<t>^60°)  for  configurations  with  blunt  noses.  In  Equation  9,  (Cp)«  =  o  is 
the  pressure  coefficient  at  «=0,  which  comes  from  Equation  3.  Figure  3  is  an 
example  of  the  application  of  Equation  9  to  a  cone  along  with  the  associated  inviscid 
surface  temperatures.  The  approximate  results  are  close  to  the  exact  cone  solution  of 
Reference  21. 
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FIGURE  3.  PERFECT-GAS  COMPARISON  OF  EXACT 
AND  APPROXIMATE  CONE  SOLUTIONS 


Figure  4  presents  the  comparison  of  the  present  methodology  for  predicting 
in  viscid  surface  temperatures  on  a  20-percent  blunt  cone  at  *  =  10  deg  and  M-r  =  15. 
These  results  are  compared  to  a  full  numerical  solution  of  the  Euler  equations 
(ZEUS)22  for  both  perfect  and  real  gases.  The  real-gas  temperatures  are  sub¬ 
stantially  lower  than  the  perfect-gas  results  and  also  agree  with  the  full  Euler 
solution  except  in  the  vicinity  of  the  overexpansion  region  past  the  blunt  tip. 
Figure  4  uses  most  of  the  theory  developed  for  the  approximate  methodology  in 
Equations  3  through  10,  along  with  the  assumptions  used  in  computing  temperature. 


2.2  AEROHEATING 

The  AP93  output  presents  boundary  layer  heating  information  in  the  form  of  a 
heat  transfer  rate,  qw;  a  heat  transfer  coefficient,  H;  and  a  recovery  temperature 
(adiabatic  wall  temperature),  Taw»  at  each  computational  point.  These  variables  are 
related  in  the  following  manner. 


aw  w 
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(spuesnoiu) 


(y)  8jnjBJ9duJ8i 
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Tw  is  the  wall  temperature.  For  high-temperature  flows,  the  heat  transfer  coefficient 
is  often  expressed  in  terms  of  enthalpies. 


At  temperatures  above  about  1500  °R,  this  enthalpy  formulation  is  the  more 
rigorously  correct  of  the  two.  Within  the  code,  all  real-gas,  high-temperature 
computations  are  performed  using  the  enthalpy  as  the  basic  variable;  but  to  maintain 
consistency  with  other  codes  in  use  at  NSWCDD,  the  results  are  converted  to  the 
temperature  form  for  output.  The  heat  transfer  is  normalized  as  shown  in 
Equations  11  and  12  because  the  coefficients  H  and  Hi  remain  fairly  constant  over  a 
wide  range  of  wall  temperatures,  even  though  the  actual  heat  transfer  rate,  qw,  may 
vary  significantly.  Thus,  since  Taw  and  haw  are  not  functions  of  wall  temperature, 
once  a  heating  computation  is  performed  for  a  given  Mach  number/altitude 
combination,  it  need  not  be  repeated  simply  because  of  changes  in  wall  conditions. 
This  weak  coupling  greatly  simplifies  the  problem  of  tracking  the  time-dependent 
thermal  response  of  a  surface  exposed  to  boundary  layer  heating.  The  aerodynamic 
solution  may  be  obtained  first  with  a  code  such  as  AP93  and  the  results  stored  in 
tabular  form  as  functions  of  Mach  number,  altitude,  and  angle  of  attack.  This 
information  can  then  be  accessed  by  an  independent  algorithm  to  compute  the  time- 
varying  heat  transfer  rates  and  the  resulting  integrated  surface  temperature  history 
along  any  given  trajectory  that  lies  within  the  limits  of  the  data  matrix. 

The  AP93  user  has  the  option,  specified  in  the  input  file,  of  using  either  perfect- 
or  real-gas  relations  in  the  computation  of  the  inviscid  solution  for  the  geometric 
configuration  and  flight  conditions  of  interest.9  In  the  former  case,  the  familiar 
analytical  thermodynamic  relations  are  employed  to  generate  the  necessary  fluid 
properties.  In  the  latter  instance,  real-gas  properties  are  determined  by  the  use  of  a 
series  of  curve  fit  relations.18.18  This  procedure  is  carried  over  into  the  routines  that 
compute  the  boundary  layer  heating.  Some  additional  work  is  required  to  generate 
the  full  complement  of  variables  needed  as  input  to  this  section  of  the  code,  and  these 
additions  are  described  in  full  in  Reference  11.  The  only  departure  from  the  use  of 
true  inviscid  surface  conditions  as  boundary  layer  edge  properties  occurs  in  the  case 
of  blunt  bodies.  The  curvature  of  the  detached  bow  shocks  associated  with  these 
configurations  creates  an  entropy  layer  near  the  body  surface.  The  inviscid  solution 
would  give  a  uniform  boundary  layer  edge  entropy  over  the  entire  body  equal  to  that 
behind  a  normal  shock  at  the  free-stream  Mach  number,  since  this  is  the  entropy 
along  the  inviscid  streamline  that  wets  the  body  surface.  In  reality,  because  of  the 
finite  thickness  of  the  boundary  layer,  the  true  edge  entropy  is  that  which  exists  at 
some  point  in  the  entropy  layer  located  at  a  distance  above  the  surface  equal  to  the 
local  boundary  layer  thickness.  This  entropy  value  is  determined  by  an  iterative 
mass  balance  technique  as  described  in  Reference  11. 
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Once  appropriate  boundary  layer  edge  conditions  are  determined,  a  series  of 
specialized  analytical  relations  is  used  to  determine  the  aerodynamic  heating  at 
various  locations.  At  the  nose  tip  stagnation  point,  a  simplified  version  of  the  Fay- 
Riddell  formula  gives23 


/dV 

q  =  0.763  Pr  06  V/pn  n  V  — -  (h  -h  ) 


aw 


w 


(13) 


The  stagnation  point  velocity  gradient,  dVe/dx,  is  determined  from  the 
Newtonian  theory,  assuming  a  spherical  nose  tip.  At  the  nose  tip,  the  flow  will 
always  be  laminar. 

If  control  surfaces  are  present,  the  viscous  heating  along  their  leading  edge 
stagnation  lines  is  determined  by  the  Beckwith  and  Gallagher  swept-cylinder 
relations24  modified  to  include  real-gas  effects. 25  For  the  laminar  case, 


qw  )  =  0.57  Pr  06\/^ 


/dve 

— —  (K  -h  )(cosA) 
dx  aw  w 


l.l 


[  14) 


where  A  is  the  leading  edge  sweep  angle  and  dVe/dx  is  the  stagnation  line  velocity 
gradient  derived  from  Newtonian  theory  assuming,  in  this  instance,  a  cylindrical 
leading  edge.  For  turbulent  flow, 


q  =  1.04  Pr 

^  W,t 


-0.6 


ip  p  > 


0.8  du  0.2 

(V  sinA)°'6(  — - )  (h  -h  ) 
on  p  dx  aw  w 


(15) 


where  Vp  is  the  flow  velocity  parallel  to  the  leading  edge  stagnation  line  and  the  (*) 
superscript  denotes  evaluation  at  a  reference  enthalpy  given  by26 


h*  =  0.5  (h  +  h  )  +  0.22  (h  -h)  (lt>) 

we  aw  e 

The  (e)  subscript  denotes  evaluation  at  the  boundary  layer  edge.  The  laminar 
or  turbulent  status  of  the  flow  is  determined  by  comparison  of  the  Reynolds  number, 
based  on  the  leading  edge  diameter,  to  user-specified  upper  and  lower  limits.  If  Ren  is 
below  the  lower  limit,  laminar  values  are  used.  If  ReD  is  above  the  upper  limit,  fully 
turbulent  flow  is  assumed.  For  intermediate  values  of  Reo,  a  linear  combination  of 
laminar  and  turbulent  values  is  output. 
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For  points  on  the  body,  the  Eckert  reference  enthalpy  flat  plate  formulation  is 
used. 27  For  laminar  flow, 


0.332  (Pr*)~a667 


and  for  the  turbulent  case, 


q 


w,t 


0.185  (Pr*)  -°-eti7 


Re* 
In - 


2.584 


(17) 


(18) 


Ni  and  Nt  are  transformation  factors  that  allow  for  the  approximation  of  three- 
dimensional  (3-D)  effects.  They  are  equal  to  three  and  two,  respectively.  The  lami¬ 
nar  or  turbulent  flow  character  is  determined  as  before  by  comparing  the  local 
Reynolds  number,  based  on  boundary  layer  running  length,  to  user-specified  upper 
and  lower  limits. 


Heating  rates  on  the  surfaces  of  wings,  fins,  or  canards  are  determined  by  using 
Equations  17  and  19,  but  in  this  case,  Nj  and  Nt  are  both  equal  to  one  because  of  the 
two-dimensional  (2-D)  nature  of  the  flow.  The  degree  of  turbulence  is  determined  in 
the  same  manner  as  for  the  body. 

Validation  of  the  new  heat  transfer  methodology  was  carried  out  in  two  steps. 
The  first  involved  comparing  results  for  stagnation  point  and  stagnation  line  flows 
from  AP93  with  those  from  the  MINIVER  code,  an  approximate  engineering  design 
tool  that  was  developed  to  model  primarily  2-D  and  axisymmetric  configurations. 
MINIVER  handles  these  specialized  situations  quite  well,  providing  a  good  check  for 
both  heat  transfer  and  inviscid  boundary  layer  edge  conditions.  The  latter 
comparisons  served  as  further  validation  of  the  real-gas  model  previously 
incorporated  into  the  AP93  code.  These  computations  covered  a  range  of  Mach 
numbers  from  7.73  to  15.  The  results  are  detailed  in  Reference  11  and  will  not  be 
repeated  here.  In  summary,  agreement  of  all  computer  variables,  including  heat 
transfer,  was  generally  good  with  variations  of  well  under  5  percent  being  typical. 
Where  larger  differences  occurred,  the  underlying  reason  was  readily  apparent,  and 
the  AP93  results  were  concluded  to  be  the  better  of  the  two.  (A  typical  case  in  point 
would  be  the  use  of  curve  fits  in  AP93  to  determine  Prandtl  number  and  viscosity, 
whereas  MINIVER  relies  on  constant  Prandtl  numbers  and  Sutherland’s  law  for 
viscosity.  At  high  temperatures,  the  curve  fits  should  be  superior.) 
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The  second  set  of  computations  in  the  validation  process  was  done  for 
spherically  blunted  cones  at  angle  of  attack.  Results  from  two  of  these  examples  are 
shown  in  Figures  5  and  6.  In  both  instances,  free-stream  pressure  and  temperature 
were  2.66  lb/ft2  and  89.971  °R,  respectively,  and  the  cone  half  angle  was  15  deg.  The 
Mach  number  was  10.6,  and  the  wall  temperature  was  540  °R.  The  results  shown  in 
Figure  5  represent  the  heating  rates  along  the  symmetry  line  of  a  cone  with  nose 
radius  of  0.375  in.  at  20-deg  angle  of  attack.  For  comparison  purposes,  analytical 
results  are  shown  from  Reference  28  that  uses  a  more  advanced  engineering  design 
code  involving  complex  streamline  tracking  methods  and  the  axisymmetric  analog 
technique  to  model  3-D  effects.  Experimental  data  from  Reference  29  are  also  shown, 
along  with  results  from  the  MINIVER  code  used  in  a  tangent  cone  mode.  AP93  and 
MINIVER  tend  to  underpredict  the  experimental  data  and  the  more  sophisticated 
computational  results  by  10  to  15  percent,  a  performance  that  is  creditable 
considering  the  simplified  solution  approach.  Figures  6  and  7  show  similar  results 
for  a  cone  with  a  1.1-in.  nose  radius  at  angles  of  attack  of  5  and  10  deg.  In  this 
instance,  comparisons  are  shown  from  AEROHEAT,  an  engineering  design  code 
similar  to  the  one  described  and  from  a  second  approximate  code3*)  that  uses  more 
sophisticated  streamline  tracking,  as  well  as  a  more  exact  boundary  layer 
formulation.  MINIVER  results  are  also  shown  along  with  experimental  data  from 
Reference  29.  AP93  does  well  beyond  about  five  nose  radii  downstream  but  tends  to 
overpredict  heating  rates  near  the  nose  tip.  This  discrepancy  is  believed  to  be  caused 
by  the  use  of  heating  expressions  that  are  based  on  flat  plate  theory  and  are  thus 
incapable  of  modeling  the  effects  of  the  rapidly  changing  boundary  layer  edge 
conditions  in  this  region.  The  AP93  results  are  slightly  better  than  those  from 
MINIVER  because  of  the  inclusion  of  variable  entropy  effects.  AP93  has  the  added 
advantage  over  MINIVER  of  being  able  to  compute  a  more  realistic  3-D  distribution 
of  heating  rates  over  the  entire  surface  at  angle  of  attack. 


2.3  SKIN  FRICTION  DRAG 

The  AP81  uses  Van  Driest  II  for  computing  skin  friction  drag  on  both  the  body 
and  fins.  The  code  assumes  a  critical  Reynolds  number  (Rec),  where  the  flow 
transitions  from  laminar  to  turbulent  of  1X106  on  the  body  and  0.5X106  on  the  fins. 
This  methodology  is  described  in  Reference  2. 

Two  changes  were  made  in  the  AP81  methodology.  These  changes  both  had  to 
do  with  the  assumptions  of  that  work.  The  AP81  methodology  was  associated  with 
typical  flight  conditions  of  vehicles  manufactured  with  a  certain  level  of  surface 
roughness  and  at  zero  angle  of  attack.  Hence,  when  comparing  the  AP81  drag  results 
to  those  of  wind  tunnel  data  (which  may  or  may  not  have  a  boundary  layer  trip)  from 
a  smooth  model,  mixed  results  were  obtained.  As  a  result,  three  additional 
alternatives  are  included  for  Rec-  The  first  is  associated  with  a  boundary  layer  trip  on 
a  wind  tunnel  model.  In  this  case,  a  small  value  for  Rec  is  assumed  so  that  the 
boundary  layer  is  turbulent  over  the  entire  model.  A  second  alternative  is  to  assume 
a  high  Rec  so  that  laminar  flow  occurs  over  the  entire  model.  A  third  alternative  is 
for  the  case  of  a  wind  tunnel  model  with  no  boundary  layer  trip.  For  that  case,  a  Rec 
of  4X106  and  2X106,  respectively,  is  assumed  for  the  body  and  wing  or  tail. 
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FIGURE  5.  WINDWARD  PLANE  HEAT  TRANSFER  RATES  FOR  0  375-IN. 
NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE  AT  «  =  20  DEG 
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FIGURE  6.  WINDWARD  PLANE  HEAT  TRANSFER  RATES  FOR  1.1  -IN. 
NOSE  RADIUS,  15-DEG  HALF- ANGLE  CONE  AT  *  =  5  DEG 
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the  options  where  «»« "»%* 

lifting  surface  ,flie'\f  ,  pceviously  discussed  at  missile  surface .  This 

from  thf  va'ueth S ^turbulent  now  occurs  over  We  ent.re^  ^  deg  to 

oc  =  30  deg  SO  that  tuny  fashion  with  angle  o 

assumption  occurs  in  a  Unea 

a  —30  deg. 


2.4  BASE  DRAG  empirical  data  for  the  body 

— r- 

alone.  Also,  an  approximate  a  ^  ^  Mach  number  based ^  tic£  ana  Space 
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Floyd  Wilcox,  coauthor  of  References  14  and  15,  was  the  chief  engineer  for  the 
tests  that  were  conducted  and  reported  in  those  references.  Eighty-nine  base 
pressure  taps  were  placed  around  a  7.2-caliber,  5-in.  diameter  body  with  a  side- 
mounted  sting.  These  taps  were  placed  every  22.5  deg  in  circumferential  location  and 
at  several  radii  from  the  body  centroid  toward  the  outer  edge.  The  configuration 
matrix  of  data  taken  is  shown  in  Table  1.  The  base  pressure  measured  at  each  of  the 
89  orifice  locations  was  then  averaged  over  its  incremental  base  area  to  get  the 
average  base  pressure  at  each  condition  of  Table  1.  Based  on  these  average  bast- 
pressure  measurements  at  each  test  condition,  changes  in  base  pressure,  and  hence, 
base  drag  because  of  a  particular  physical  model  change  or  flight  condition  change 
could  be  readily  computed  by  simply  substracting  the  two  data  points. 

TABLE  1.  CONFIGURATION  INDEX 


Config 

Fins 

Off 

t/c 

0.05 

0.10 

015 

0.5 

x/c 

1.0  1.5  2.0 

25 

0 

5 

s 

10 

15  20 

at 

<M*=2.0) 

(M,  2  2.51 

1 

X 

Sweep 

Sweep 

2 

X 

X 

X 

0,5,10 

0 

3 

X 

X 

X 

0,5,10 

0 

4 

X 

X 

X 

0,5,10 

0 

5 

X 

X 

X 

0,5,10 

0 

6 

X 

X 

X 

0,5,10 

0 

7 

X 

X 

X 

0.5,10 

0 

8 

X 

X 

X 

0,5,10 

0 

9 

X 

X 

X 

0,5.10 

0 

10 

X 

X 

X 

0,5,10 

0 

11 

X 

X 

X 

0,5,10 

0 

12 

X 

X 

X 

0.5,10 

0 

13 

X 

X 

X 

0,5,10 

0 

14 

X 

X 

X 

0,5.10 

0 

15 

X 

X 

X 

0,5,10 

0 

16 

X 

X 

X 

0,5,10 

No  data 

Using  the  process  described,  along  with  a  wind  tunnel  database  not  available 
when  AP81  was  developed, 31  a  new  empirical  estimate  of  Kase  pressure  coefficient 
(CP|!)  was  derived.  This  new  estimate  is  shown  in  Figure  8  and  compared  to  the  AP81 
value  of  CPli.  The  two  curves  are  similar,  with  the  AP93  slightly  higher  than  AP81 
for  M-r  <15  and  slightly  lower  than  AP81  for  M*-2:3.0.  Body-alone  angle-of-attack 
effects  on  base  pressure  are  then  estimated  by 


NF. 


-0 


1  -FO.OtFj 


(19) 
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Here,  (CPl!)s,f  _0  comes  from  Figure  8  and  Ft,  the  increase  due  to  angle  of  attack  from 
Figure  9a  or  9b,  depending  on  the  Mach  number.  Boattail  and  power-on  effects  on 
base  drag  are  estimated  as  present  in  AP81. 


FIGURE  8.  MEAN  BODY-ALONE  BASE  PRESSURE 
COEFFICIENT  USED  IN  AP81  AND  AP93 


At  this  point,  it  is  worth  noting  that,  while  the  databases  of  References  14,  15, 
and  31  helped  to  improve  the  estimate  of  base  pressure  as  a  function  of  Mach  number 
and  angle  of  attack  for  the  body  alone,  additional  data  are  still  needed  for  «  >  15  deg 
at  all  Mach  numbers.  This  need  is  indicated  by  the  dotted  lines  in  Figures  9a  and  9b, 
which  are  extrapolations  from  data  available  for  «<15  deg  and  engineering 
judgment.  This  same  statement  will  also  be  even  more  true  for  fin  effects  due  to 
control  deflection  and  angle  of  attack,  as  will  be  discussed  in  the  following  para¬ 
graphs. 

The  total  body  base  pressure  coefficient  for  fins  located  flush  with  the  base  is 

(c  ) 

^B' *.fi,t/c,x/c  =  0 

where  (Cpn)N(,_0,  F2,  and  F3  come  from  the  AP93  curve  of  Figures  8,  10,  and  11, 
respectively. 


1+0. OIF 


•  <CJ, 


+  0.01  F.j't/d 
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FIGURE  9.  PERCENT  INCREASE  IN  BODY-ALONE  BASE  PRESSURE 
COEFFICIENT  DUE  TO  ANGLE  OF  ATTACK 
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In  Figure  10,  no  data  were  taken  for  M*,<2  in  the  test  of  References  14  and  15, 
and  none  could  be  found  in  the  literature.  Hence,  the  data  for  M-r>  =  2  are  assumed  to 
apply  for  Ma><2  as  well.  While  this  is  a  big  assumption,  it  is  believed  to  be  better 
than  neglecting  the  base  pressure  effect  due  to  control  deflection  and  angle  of  attack, 
which  other  engineering  aerodynamics  codes  do.  It  is  also  worth  noting  that  Figure 
11  indicates  what  is  intuitively  obvious;  i.e.,  for  small  control  deflections  and  angles 
of  attack,  fin  thickness  effects  are  important  in  base  pressure  estimation,  whereas  for 
large  values  of  «  and  8,  the  additional  change  in  Cp(1  due  to  fin  thickness  is  minimal. 


The  final  parameter  to  define  the  effect  on  base  pressure  is  fin  location  relative 
to  the  body  base.  This  is  done  through  Equation  21,  where 


(CJ 


B  ■  a.fi.t/c.x'c 


+  0.01 


b7nf.* 


(acp  ) 

rR  /  ,  K 


(21) 


B '  *.fi,t/c,xe 


Here  (Cpu)NKa  is  the  body-alone  base  pressure  coefficient  at  a  given  angle  of  attack 
given  by  Equation  19  and  (ACpu)^  8t/cx/c  is  the  total  change  due  to  the  presence  of  fins 
at  a  given  «,  8,  t/c,  and  x/c.  An  example  of  (ACp„),fit/cx/c  is  given  in  Figure  12  for 
M*,  =  2.0  and  l<* +  81  =  10  deg.  Reference  14  shows  other  curves  for  this  parameter. 
Figure  12  shows  that  the  change  in  base  pressure  due  to  all  variables  present  varies 
from  that  at  x/c  =  0,  where  the  fins  dominate  to  that  of  the  body  alone  where  the  fins 
have  no  effect  (x/c  =»  2.5). 


FIGURE  10.  PERCENT  INCREASE  IN  BASE  PRESSURE  COEFFICIENT 
DUE  TO  COMBINED  EFFECTS  OF  ANGLE  OF  ATTACK 
ANDCONTROL  DEFLECTION  (t/c  =  0> 
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DATA  SOURCES 


M, 


FIGURE  11.  ADDITIONAL  PERCENT  INCREASE  IN  BASE  PRESSURE 
COEFFICIENT  DUE  TO  FIN  THICKNESS 
AT  VARIOUS  VALUES  OF  I*  +  81 


X/C 


FIGURE  12.  PERCENT  INCREASE  IN  BASE  PRESSURE  COEFFICIENT 
DUE  TO  FIN  LOCATION  I*  +  81  =  1 0  DEG,  IVU  =  2.0 
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2.5  BODY-ALONE  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 
The  normal-force  coefficient  of  the  body  alone  is  estimated  by 


CN  =  CN, 


+  c 


N 


NL 


(22) 


where  Cn,  is  the  linear  term  and  Cnm  the  nonlinear  term.  The  linear  term  is 
predicted  in  AP81  by  either  SOSET,  Second-Order  Van  Dyke  combined  with  MNT,  or 
empirical  depending  on  the  Mach  number  range  (see  References  2,  7,  or  9  for  details). 
The  nonlinear  term  is  estimated  by  the  Ailen-Perkins  viscous  crossflow  theory  34  No 
changes  were  made  in  the  linear  term  of  Equation  22  in  AP93  from  AP81.  Three 
change^  in  the  nonlinear  term  of  Equation  22  were  made  for  the  AP93. 


The  nonlinear  term  of  Equation  22  according  to  Reference  34  is 


A 

C.,  =  rjC,  sin2«  — - 

NNL  dc  A 


ref 


(23) 


The  first  change  from  AP81  is  in  the  value  of  q.  AP81  used  an  incompressible 
value  of  q  with  no  account  of  compressibility  effects.  Reference  35  clearly  showed  a 
compressibility  effect.  This  effect  is  shown  in  Figure  13a  along  with  the  line  drawn  to 
represent  the  data.  This  line  is  defined  as 


n  = 


1  — q 
*o 

1.8 


M  M  +  q  for  ^  1.8 

N  lo  N 


(24) 


n  =  l 


for  M  N  >  1.8 


where  q0  is  the  incompressible  value  of  q  (Mn  =  0)  given  in  Reference  2  and  used  in 
AP81. 


The  second  change  is  in  the  value  of  the  crossflow  drag  coefficient  used.  This 
value  was  changed  to  allow  the  effect  of  transition  on  the  body  surface  to  impact  the 
value  chosen.  This  impacts  the  value  of  Cdc  for  Mn  values  of  0.5  and  less.  Also,  the 
value  of  Cdc  is  slightly  lower  for  0.6  sMn  ^2.2  than  that  used  in  AP81.  This  is  based 
on  the  large  NASA  Tri-Service  Database. 36  The  new  value  of  Cdc  used  in  AP93  is 
given  in  Figure  13b.  If  the  flow  on  the  body  is  a  combination  of  laminar  and 
turbulent  (which  is  the  case  for  most  conditions),  a  value  somewhere  in  between  the 
two  values  on  the  Figure  13b  curve  for  Mn  ^0.5  will  be  computed.  If  XL  defines  the 
length  of  laminar  flow  on  the  body  and  Xx  is  the  total  length,  then  for  Mn^0.5, 


C,  =  1.2  - 

dc 


(25) 


Thus,  if  Xl  =  0  so  the  flow  over  the  body  is  fully  turbulent,  a  value  of  Cdc  — 12  will  be 
computed,  whereas  a  value  of  0.4  will  be  picked  if  the  flow  is  fully  laminar. 
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(MN  =  M«SINa) 

FIGURE  13a.  COMPRESSIBILITY  EFFECTS  ON  CROSSFLOW 
DRAG  PROPORTIONALITY  FACTOR 


FIGURE  13b.  CROSSFLOW  DRAG  COEFFICIENT  FOR 
AN  OGIVF -CYLINDER  CONFIGURATION 


FIGURE  13c.  CENTER  OF-PRESSURE  SHIFT  IN  BODY-ALONE 
NORMAL  FORCE  FOR  « > 10 DEG 
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The  third  change  made  in  AP93  was  in  the  center-of-pressure  location.  AP81 
used  a  weighted  average  of  the  normal  force  center  of  pressure  of  the  linear  term  and 
nonlinear  term,  where  the  nonlinear  term  Xcp  was  at  the  centroid  of  the  planform 
area  in  the  crossflow  plane  and  the  Xcp  of  the  linear  term  was  computed  theoretically 
or  empirically.  Both  these  values  were  held  constant  as  angle  of  attack  increased,  the 
only  change  being  from  the  changing  values  of  the  normal-force  terms  of  Equation  22. 
In  numerical  experiments  using  the  database  of  Reference  36,  it  was  found  that  the 
assumption  of  a  constant  value  of  center  of  pressure  with  angle  of  attack  was  not 
completely  correct.  It  is  suspected  that  as  angle  of  attack  increases,  the  center  of 
pressure  of  the  linear  term  of  Equation  22  changes  and  can  no  longer  be  assumed  to 
be  constant.  An  empirical  way  to  represent  this  change  with  Mach  number  is  given 
in  Figure  13c  This  change  is  effective  for  *  >  10  deg.  Between  *  =  0  and  10  deg,  the 
correction  is  implemented  in  a  linear  fashion  between  zero  at  *  =  0  to  its  full  value  at 
«  =  10  deg. 

Figure  14  is  an  example  of  the  normal-force  and  center-of-pressure  comparisons 
of  the  AP81,  AP93,  and  experimental  data.  The  data  are  from  Reference  36,  which  is 
for  a  12.33-caliber  tangent-ogive  cylinder  configuration  with  a  3.0-caliber  nose.  The 
improvements  made  in  AP93  give  significantly  better  results  on  both  CN-  and  Xcp  as  a 
function  of  angle  of  attack. 


2.6  WING-ALONE  LINEAR  NORMAL  FORCE  AND  CENTER 
OF  PRESSURE  FOR  LOW-ASPECT  RATIO  WINGS 

The  AP81  gives  reasonably  accurate  wing-alone  linear  normal-force  coefficients 
if  the  aspect  ratio  is  one  or  greater.  However,  as  the  aspect  ratio  gets  small,  similar  to 
what  would  be  obtained  on  dorsal-like  lifting  surfaces,  the  linear  theory  gives  values 
of  (CnJ^q,  which  are  too  high.  This  is  because  the  AP81  neglects  Mach  lines  that 
intersect  opposite  wing  tips.  This  assumption  becomes  increasingly  invalid  as  the 
aspect  ratio  gets  small.  To  remedy  this,  the  Evvard^7  and  Krasilshchikova38 
methods  for  low-aspect  ratio  wings  were  implemented  in  the  AP93.  The  (C.\x),  _0  is 
still  slightly  higher  than  experimental  data;  however,  it  is  believed  this  is  due  to 
neglecting  thickness  effects.  For  the  details  of  this  effort,  interested  readers  are 
referred  to  Reference  12. 

A  second  change  was  also  made  in  the  linear  part  of  the  wing-alone  lift.  This 
had  to  do  with  center-of-pressure  prediction.  When  the  linear  theory  or  lifting 
surface  theory  said  the  center  of  pressure  was  further  forward  on  the  wing  than 
15  percent  {or  10  percent  of  the  chord  subsonically),  the  AP93  nuts  a  limit  at  these 
values.  This  limit,  in  effect,  brings  into  the  center-of-pressure  considerations  the  real 
viscous  effects,  which  are  neglected  in  the  linearized  inviscid  theories. 
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Angle  of  Attack  (deg) 


FIGURE  14.  BODY-ALONE  NORMAL-FORCE  COEFFICIENT 
AND  CENTER  OF  PRESSURE  (M»  =  3.5) 
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2.7  WING-ALONE  NONLINEAR  NORMAL  FORCE 
AND  CENTER  OF  PRESSURE 


One  of  the  major  reasons  the  AP8I  gave  poor  results  at  <*>10  deg  for  many 
missile  configurations  was  failure  to  include  nonlinearities  in  the  wing  lift.  Using 
the  databases  of  References  39  and  40,  a  semiempirical  method  was  developed  for  the 
nonlinear  wing-alone  normal-force  term!2,l3  analogous  to  the  body-alone 
Equations  22  and  23;  i.e.,  the  nonlinear  term  of  wing-alone  lift  can  be  defined  as 


=  f(MN,AR,A) 


(26) 


Here,  f(Mu,  AR,  A)  is  analogous  to  the  q  Cdc  of  the  body  alone  in  Equation  23.  Since 
the  total  wing-alone  normal  force  is  known  for  a  given  AR,  M*,  A,  and  «  from 
References  39  and  40,  the  linear  value  of  lift  is  known  from  the  3-D  thin-wing  theory 
or  lifting  surface  theory  from  AP81  (modified  according  to  Section  2.6);  the  nonlinear 
normal  force  of  the  wing  alone  is 

CN  (Mn,  AR,  A)  =  CN  <Mn,  AR,  A)  -  CN  (M  ,  AR,  A)  (27) 

NL  L 

Using  the  data  of  References  39  and  40,  Equation  27  values  were  generated  and 
a  parameter  ki  defined  as 


k 


L 


C  (Mn,AR,A) 

NL 

sin2* 


(28) 


was  generated.  Tables  of  ki  for  both  high  and  low  Mach  numbers  are  given  in 
Tables  2  and  3.  The  total  wing-alone  normal  force  in  AP93  is  therefore 


+  kj  sin2* 


(29) 


The  second  term  of  Equation  29  was  neglected  in  AP81. 


The  center  of  pressure  of  the  wing-alone  lift  was  assumed  to  vary  quadratically 
between  its  linear  theory  value  at «  =  0  to  the  centroid  of  the  planform  area  (adjusted 
for  thickness  effects)  at  *  =  60  deg. 


Defining  the  center  of  pressure  of  the  wing-alone  linear  term  as  A  and  the 
center  of  pressure  of  the  nonlinear  term  as  B  (both  in  percent  of  mean  geometric 
chord),  then  the  center  of  pressure  of  the  wing  lift  is 


(X  ) 


cp  w 


=  A  + 


—  I*  I (B~ A|  + 
36  w 


1 

5400 


«2  |A-B| 

W 


(30) 
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TABLE  2.  VALUES  OF  kt  FOR  LOW  MACH  NUMBERS 


a/m* 

0.0 

0.5 

1.0 

0.0 

1.55 

1.57 

1.60 

0.5 

2.84 

2.90 

2.82 

1.0 

2.37 

2.45 

2.43 

2.5 

3.0 

3.5 

4.0 

4.5 

1.25 

0.92 

0.56 

0.29 

0.16 

1.00 

0.80 

0.64 

0.47 

0.33 

1.05 

0.90 

0.75 

0.61 

0.48 

0.5 

1.0 

1.5 

1.48 

1.46 

0.99 

2.45 

1.85 

0.70 

1.22 

1.10 

0.50 

AR  =  1.0;  M^  <  3.5 


2.0 


2.5 

3.0 

3.5 

4.0 

0.22 

0.12 

0.09 

0.09 

0.19 

0.20 

0.26 

0.36 

0.50 

0.65 

0.78 

0.88 

AR  >  2.0;  Mx  <  3.5 


0.5 

1.0 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

-1.84 

-1.95 

-1.50 

-0.20 

0.00 

0.10 

0.20 

mvm 

0.30 

-1.84 

-1.95 

-1.50 

-0.20 

0.30 

0.41 

0.60 

pus 

0.80 

-1.47 

-1.35 

0.70 

0.20 

0.60 

0.83 

0.98 

1.09 

1.15 

FIGURE  3.  VALUES  OF  kt  FOR  HIGHER  MACH  NUMBERS 
AR  <  0.5;  M,  >  4.0 


A / 

Mocsin® 

0.0 

0.5 

1.0 

1.5 

0.0 

-1.60 

-0.98 

0.23 

0.55 

0.5 

-0.87 

-0.24 

0.33 

0.60 

1.0 

-0.31 

0.09 

0.46 

0.68 

m 

Be  a 


2.5 

3.0 

3.5 

0.82 

0.89 

0.92 

0.82 

0.89 

0.92 

0.87 

0.91 

0.93 

.5  5.0  5.5  6.0 


0.95  0.95  0.95  0.95 


0.95  0.95  0.95 


AR  =  1.0;  Mx  >  3.5 


A I 

Musin'* 

0.0 

0.5 

1.0 

H 

2.0 

2.5 

3.0 

3.5 

4.0 

m 

5.0 

5.5 

6.0 

0.0 

-0.39 

-0.39 

-0.29 

0.06 

0.29 

0.48 

0.60 

0.69 

0.75 

0.81 

0.86 

0.91 

0.5 

pjjjj 

0.17 

0.29 

0.46 

0.63 

0.76 

0.85 

0.90 

0.93 

0.95 

0.95 

mm 

1.0 

EEEI 

0.50 

0.86 

0.93 

0.94 

0.95 

0.95 

0.95 

0.95 

ififfil 

0.95 

0.95 

0.95  | 

AR  >  2.0;  M*  >  3.5 


K) 

M^sin* 

0.0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

H 

5.0 

5.5 

6.0 

0.0 

-0.25 

-0.05 

0.20 

0.50 

0.80 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.5 

0.02 

0.29 

0.80 

0.98 

EE9 

0.97 

0.97 

0.95 

0.95 

0.95 

EES 

1.0 

0.66 

1.02 

1.15 

1.18 

1.09 

1.02 

0.95 

0.95 

0.95 

EE3 

Em 
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«w  in  Equation  30  is  the  total  angle  of  attack  in  degrees  on  the  wing.  Figure  15  gives 
an  example  of  the  AP93  methodology  compared  to  AP81  and  experimental  data.  This 
particular  case  shows  significant  improvement  in  the  wing-alone  normal  force  of  the 
AP93  versus  AP81  when  compared  to  the  experiment.  However,  no  improvement  in 
center  of  pressure  is  obtained  because  A  =  0  and  the  centroid  of  the  planform  area  are 
the  same  as  experimental  data  suggest. 


FIGURE  15.  WING-ALONE  NORMAL-FORCE  COEFFICIENT  AND 
CENTER  OF  PRESSURE  <  AR -  0.5,  A  =  0.0,  =  1.6) 
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2.8  WING-BODY  AND  BODY-WING  NONLINEAR  INTERFERENCE 
FACTORS  DUE  TO  ANGLE  OF  ATTACK 

The  total  configuration  normal-force  coefficient  at  a  given  angle  of  attack, 
control  deflection,  and  Mach  number  is41 


The  first  term  in  Equation  31  is  the  normal  force  of  the  body  alone  including  the 
linear  and  nonlinear  components;  the  second  term  is  the  contribution  of  the  wing  (or 
canard)  including  interference  effects  and  control  deflection;  the  third  term  is  the 
contribution  of  the  tail  including  interference  effects  and  control  deflection;  and  the 
last  term  is  the  negative  downwash  effect  on  the  tail  due  to  wing  shed  or  body  shed 
vortices.  The  K’s  represent  the  interference  of  the  configuration  with  respect  to  angle 
of  attack,  and  the  k’s  represent  the  interference  with  respect  to  control  deflection. 
Each  of  these  interference  factors  is  estimated  in  the  AP81  by  slender  body  or  linear 
theory  .3  As  such,  they  are  independent  of  angle  of  attack.  Estimating  the  change  in 
these  parameters  with  angle  of  attack  is  desired.  To  do  this,  Reference  12  made  use  of 
the  large  database  of  References  36, 39,  and  40. 

The  work  of  Reference  12  found  that  the  wing-body  interference  factor  had  the 
qualitative  behavior  as  shown  in  Figure  16.  At  low  angles  of  attack,  slender-body 
theory  appeared  to  be  a  good  estimate  of  KW(B).  This  estimate  was  adjusted  slightly 
for  Moo ~1.5  by  an  amount  AKwib>-  At  some  angle  of  attack  defined  as  <xc,  KW(B) 
seemed  to  decrease  in  a  nearly  linear  fashion.  The  rate  of  this  decrease  was  a 
function  of  Mach  number — the  higher  the  Mach  number,  the  larger  the  rate  of 
decrease.  At  some  point  defined  as  <*q,  the  Kw<bi  appeared  to  reach  a  minimum  and 
remain  about  constant.  As  a  result  of  this  analysis,  a  mathematical  model  was 
derived  to  define  KW(Bi  in  terms  of  its  slender-body  theory  value  [KW(B)lsB  and  an 
empirical  correction  derived  from  the  databases  of  References  36,  39,  and  40.  This 
model  given  in  Figure  16  is 


^W(B) 


*SviB>  SB  +  AKW< 
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■O'  0.5  1 
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W(B) 
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i  =  0 


dKWtB>  1  r/s 

4.  -  (a  — a  )>  -  for* 

d«  c  1  0.5  c  ® 
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GENERAL  EQUATION: 


W(B) 


WcB) 


Kw»'i»  +1AK*»>L-.o5 

akwJ  I !Lj*Jhs£L)<n?™lmm 

w<m  io-oO.5  V  jo  /  0.5  “  * 


acaca- 


KW.»,  *  +||AKW(»>|..4  +  ~^!1(0D-“,>|(^)fi'o>«D 

FIGURE  16.  QUALITATIVE  BEHAVIOR  OF  WING-BODY 
INTERFERENCE  FACTOR  AS  A  FUNCTION 
OF  ANGLE  OF  ATTACK 


The  empirical  corrections  to  Kw<b)  are  also  in  a  form  that  can  be  defined  mathe¬ 
matically  as  opposed  to  a  table  lookup  procedure.  These  equations  for 


AKw 


(B) 


dKW(B> 

„ ,  - ,  «  ,  r  nd  » 

*  =  0  d<x  c  D 


are 


AK 


W(B) 
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I  AKW(B) 

=  0.22 

-*  =  0 

( AKW(B>  j 

=  -0.44 

*  =0 

M  -  1.5 

an 
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=  0 


for  M  S  1.0 


for  1.0  <  M  ^1.5 


for  M  >1.5 
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d 

/  d« 

W(B) 

d<* 


-(0.00283M  +0.025) 

QO 


(34) 


QC 

C 

M<2.0 

x  =12.5-  1.06  M  -  2.59  M2  forAR<0.5 

C  or  or 

oc  =12.5- 6.25  M  for  AR  =  1.0 

c  * 

=  4.5  +  2.25  Mx  -  2.25  for  AR>2.0 

Moo  >2.0 
«c  =  0 


D 


<*  =  33.3  -  8.19  M  +  0.82  M 

D  JO  a 


for  A  =  0 


«_  =  25.3  -  6.62  M  +0.66M" 

D  a-  a 


t,r  A  =  1.0 


“d  ^D^=t.0  +  ^ 


'“d^O  **D*A  =  1.0 


for0<A<1.0 


(35) 


(36) 


The  semiempirical  model  for  Kg,w,  was  also  defined  in  terms  of  its  slender  body 
or  linear  theory  value,  plus  a  correction  due  to  nonlinearities  associated  with  angle  of 
attack.  The  mathematical  model  for  Kgiwi  was  defined  as*2 


iW> 
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r/s  f 
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B(W) 
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I.T 

K 


B<W> 
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Unfortunately,  a  mathematical  model  for  (AKbiw.<L=o  and  d|KB,w>J/da  were  difficult 
to  define  because  of  the  variability  of  the  constants  as  a  function  of  the  parameters  of 
interest.  As  a  result,  a  three-parameter  table  lookup  for  these  two  parameters  is  used 
in  AP93  based  on  the  data  in  Table  4.  The  parameters  in  the  table  lookup  include 
Moo,  A,  and  AR.  Linear  interpolation  is  used. 

TABLE  4.  DATA  FOR  BODY-WING  NONLINEAR  SEMIEMFIRCAL 
INTERFERENCE  MODEL 


Aspect 

Ratio 

Taper 

Ratio 

£0.6 

0.8 

Data  for  (  AKg,-^  ]  Q=o 

Mach  Number 

1.2  1.5  2.0  25 

3  0 

3.5 

£4.5 

<0.25 

0,0.5, 

1.0 

-0.1 

-0.1 

0.5 

0.6 

0.7 

0  8 

0.7 

0.5 

0.3 

0.5 

0.5 

-0.28 

-0.1 

0.13 

0.11 

0.05 

-0.02 

-0.06 

0 

0 

1.0 

0.5 

-0.26 

-0.2 

0.15 

0.21 

0.15 

0 

0 

0 

0 

£2.0 

0.5 

-0.13 

-0.04 

0.12 

0.43 

-0.16 

0 

0.37 

-0  08 

-0  16 

0.5 

0 

•0.3 

-0.06 

0.26 

0.28 

0.17 

0.12 

0  14 

0 

0 

£2.0 

0 

-0.2 

-0.1 

0.12 

0.52 

0.12 

0.15 

0.22 

006 

-0.22 

0.5 

1.0 

-0.16 

0.08 

0.26 

0.14 

-0.12 

0 

-0.05 

-0.10 

0 

£2.0 

1.0 

-0.2 

-0.1 

0.12 

0.45 

-0.02 

Oil 

0  28 

-0.17 

-0.3 

Aspect 

Ratio 

Taper 

Ratio 

£0.6 

0.8 

Data  for 

1.2 

d[  Kb;w>  )  1  da 

Mach  Number 

1.5  2.0  2.5 

3.0 

3.5 

£4  5 

£0.25 

0,0.5, 

1.0 

0.018 

0.013 

-0.010 

-0.023 

-0.013 

-0.022 

-0.031 

-0  025 

-0.031 

0.5 

0.5 

0.019 

0.010 

-0.008 

-0.010 

-0.013 

-0.013 

-0.013 

-0  012 

-0  012 

1.0 

0.5 

0.013 

0.010 

-0.007 

-0.013 

-0.020 

-0.017 

-0.012 

-0  012 

-0.012 

£2.0 

0.5 

0.010 

0.011 

0 

-0.013 

-0.010 

-0.017 

-0  040 

-0.012 

-0  012 

0.5 

0 

0.033 

0.022 

0 

-0.007 

-0010 

-0.008 

-0.014 

-0.012 

-0.012 

£2.0 

0 

0.010 

0.010 

-0.007 

-0.020 

-0.011 

-0.020 

-0.023 

-0.012 

-0  012 

0.5 

1.0 

0.019 

0 

-0.019 

-0.010 

-0.007 

-0.013 

-0.014 

-0.012 

-0  012 

£2.0 

1.0 

0.010 

0.01 

-0.007 

-0.017 

0 

•0.017 

-0.026 

-0  012 

-0.012 

In  Equations  32  and  36,  the  factor 


r/s 

0.5 

appears.  This  is  because  the  database  of  Reference  36  is  based  on  r/s  =  0.5,  and 
Reference  41  indicates  the  aerodynamics  vary  linearly  with  r/s.  This  assumption  is 
inherent  in  the  semiempirical  models  for  Kw/b>  and  Kb,w.- 
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In  examining  cases  where  r/s  is  small,  it  was  found  that  at  high  angles  of  attack, 
the  wing-alone  solution  was  not  recovered  properly  through  the  Equations  32  and  37 
process.  To  remedy  this  situation,  the  AP93  nonlinear  interference  factors  were 
blended  into  those  predicted  by  slender-body  or  linear  theory  as  r's  became  small. 
The  specific  equations  used  to  do  this  are 

For  r/s  ^0.25 


^WiBi  ^WlB> 
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For  0.05  <  r/s  <0.25 
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SBT  ^ 
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BiW) 

SBT 
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j  (r/s-  0.05i/0.2 


(38b) 


For  r/s  <0.05 
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;  K 


SBT 


BiW) 


K 


Bt  Wi 


SBT 

I.T 


(38c) 


In  essence,  the  model  represented  by  Equations  38a  through  38c  uses  the 
nonlinear  interference  factors  for  r/s  values  greater  than  0.25;  they  use  a  blend  of 
slender-body  or  linear  theory  and  the  nonlinear  values  of  interference  factors  for  r/s 
values  between  0.05  and  0.25.  They  also  use  the  slender-body  or  linear  theory  values 
for  r/s  values  less  than  0.05.  Hence,  when  the  body  vanishes  (r/s  =  0),  the  wing-alone 
solution  will  be  automatically  recovered  in  a  smoother  and  more  accurate  way. 


Figure  17  is  an  example  of  the  normal  force  on  the  wing  in  the  presence  of  the 
body  and  the  normal  force  on  the  body  in  the  presence  of  the  wing  using  the  AP93 
theory,  the  AP81  theory,  and  compared  to  experimental  data.  Note  that 


CN  =  CN  KW(B) 
Wi  Bi  W  D 


Cv  =  Cv  K 

NBIW>  nw  b,Wi 


(39) 
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Mac* -1.2 


0  S  10  15  20  25  30 

Angle  of  Attack  (deg) 


Mach-12 


0  5  10  IS  20  23  30 

Angle  of  Attack  (deg) 


FIGURE  17.  WING-BODY  AND  BODY-WING  INTERFERENCE  AS 
A  FUNCTION  OF  «  ( AR  =  2.0,  A  =  0,  M*  =  1.2) 


Hence,  Figure  17  is  actually  a  representation  of  the  normal-force  coefficient  on  the 
wing  and  additional  normal  force  on  the  body  due  to  the  wing.  Thus,  Equation  39  is  a 
representation  of  the  accuracy  of  not  only  Kwibi  and  Kb<w),  but  Cnw  in  conjunction 
with  the  interference  factors.  This  is  a  more  true  indication  of  the  accuracy  of  the 
code,  because  there  are  actually  two  of  the  component  force  terms  that  make  up 
Equation  39.  As  seen  in  Figure  17,  the  AP93  methodology  is  clearly  superior  to  the 
AP81  theory  as  angle  of  attack  increases. 


2-30 


NSWCDD/TR-93/91 


The  center  of  pressure  of  the  new  value  of  normal  force  of  the  wing  in  the 
presence  of  the  body  estimated  by  Equation  32  is  assumed  to  remain  at  the  values  of 
the  wing-alone  solution  of  AP93  given  by  Equation  30.  The  center  of  pressure  of  the 
additional  lift  on  the  body  due  to  the  presence  of  the  wing  is  estimated  the  same  as  in 
AP81,  which  is  either  slender-body  or  linearized  theory.  These  values  are  modified 
for  short  afterbodies  according  to  Reference  3. 

In  exercising  the  AP93  on  missile  configurations  in  the  transonic  speed  regime 
(0.6<M<2.0),  it  was  found  that  some  of  the  nonlinear  lift  associated  with  small 
aspect  ratio  fins  (AR<  1.4)  was  lost  due  to  shock-wave  formation.  A  certain  amount 
of  linear  lift  loss  was  accounted  for  from  an  empirical  approach  in  the  AP81.  This 
appeared  to  be  satisfactory  for  the  larger  aspect  ratio  fins,  where  the  nonlinear 
normal-force  term  with  angle  of  attack  was  negative.  However,  when  the  fins  have  a 
positive  nonlinear  normal  force  due  to  angle  of  attack,  some  of  this  force  appears  to  be 
lost  with  shock  waves.  This  loss  was  estimated  empirically  as  a  function  of  Mach 
number  and  angle  of  attack  for  a  wing  that  had  an  area-to-body  reference  area  of 
about  one.  These  data  for  AC.\  losses  due  to  compressibility  effects  are  given  in 
Table  5.  A  two  parameter  linear  interpolation  is  made  from  Table  5  for  a  given  Mr 
and  «  to  compute  ACn.  ACn  is  further  degraded  for  taper  ratio  for  values  of  A <0.5. 
The  specific  equations  for  AC*  are 

/  \ 

ACV  =  -  (  ACV.  -  for  A  >0.5 

NB.W.  1  N-/Aref 


TABLE  5.  LOSS  OF  WING  NONLINEAR  NORMAL  FORCE  DUE  TO 
SHOCK-WAVE  EFFECTS  IN  TRANSONIC  FLOW 
|*  +8|  DEG 


Mac 

0 

5 

10 

15 

20 

25 

30 

35 

IV 

o 

<0.4 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.6 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0220 

-0.2060 

-0.6890 

-0  9500 

-1.300 

0.8 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0531 

-0.2200 

-0.7100 

1  010 

-1.400 

1.2 

0.0000 

0.0000 

-0.0093 

-0.0293 

-0.1651 

-0.4167 

-0.7629 

-1.070 

-1  500 

1.5 

0.0000 

0.0000 

-0.0653 

-0.1111 

-0.1556 

-0.4444 

-0.7000 

-1.070 

-1.500 

2.0 

0.0000 

0.0000 

-0.0076 

-0.0376 

-0.1502 

-0.1142 

-0.0951 

-0.0700 

-0.0500 

>2.5 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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2.9  NONLINEAR  WING-BODY  INTERFERENCE  FACTOR 
DUE  TO  CONTROL  DEFLECTION 

The  nonlinear  methodology  for  computing  the  wing-body  interference  due  to 
control  deflection  has  not  been  presented  before.  As  a  result,  a  more  detailed  dis¬ 
cussion  of  this  effort  will  be  presented  in  contrast  to  the  other  new  methods,  which 
were  only  briefly  summarized. 

Initially,  it  was  planned  to  use  slender-body  theory  for  the  interference  factors 
kW(B)  and  kB(wt,  as  currently  done  in  AP81.  This  plan  was  based  on  Reference  12 
results  comparing  computations  (using  Equation  31  where  all  the  nonlinearities  are 
included)  with  experimental  data  at  8  =  0  for  both  body -tail  and  body -wing-tail  or 
body-dorsal-tail  configurations.  These  comparisons  were  good  and  seemed  to  indicate 
this  new  technology  to  be  superior  to  existing  engineering  approaches.  However, 
when  results  were  examined  for  configurations  that  had  control  deflections  on  either 
the  aft  or  forward  lifting  surface,  they  were  found  to  not  be  as  good  as  desired.  This 
led  to  the  conclusion  that  nonlinear  interference  factors  due  to  control  deflection  were 
also  required  to  improve  the  performance  of  AP93  when  compared  to  experimental 
data. 

To  address  this  problem,  the  database  of  Reference  36  was  once  again  employed. 
An  attempt  was  made  to  estimate  kw,B>  based  on  the  component  measurements  of  this 
database.  During  these  measurements,  only  one  fin  was  deflected  at  a  given  set  of 
conditions  (*,  M®,  <|>,  fin  geometry).  Referring  to  Equation  31,  the  normal-force 
coefficient  component  resulting  from  a  control  deflection  of  a  lifting  surface  is 


AC 


N 


W(B) 


<Bi 


(CN 


W 


W 


(41) 


The  additional  normal-force  coefficient  on  the  body  due  to  control  deflection  is 


AC‘nb<w)  kB'W>(cNJw5w  (42) 

In  principal,  because  the  database  of  Reference  36  was  taken  from  wing  and 
body  independently  where  they  were  close  enough  to  each  other  to  get  the 
interference  effects,  kW{B)  and  kB<w)  can  be  determined.  This  was  attempted  with 
mixed  success.  Values  of  single-fin  control  deflection  normal-force  coefficients  were 
obtained  for  several  conditions.  The  ACNwiK  of  Equation  41  could  be  readily 
determined  by  multiplying  the  values  in  the  wind  tunnel  data  by  two  (since  only  one 
fin  was  deflected)  and  subtracting  these  values  for  the  AC.\Wi(U  at  8  =  0. 

It  is  suspected  that  there  are  several  reasons  for  the  mixed  success  of  this 
approach.  First  of  all,  the  area  of  any  one  fin  that  was  deflected  varied  in  physical 
size  from  0.0312  to  0.0078  ft2.  This  compared  to  a  body  cross-section  area  of  0.0490  ft2 
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and  a  body  planform  area  in  the  crossflow  plane  of  approximately  0.7  ft2.  In  other 
words,  at  zero  angle  of  attack,  where  body  lift  and  wing  lift,  due  to  angle  of  attack  are 
zero,  a  reasonable  value  for  the  normal-force  coefficient  resulting  from  control 
deflection  can  be  obtained.  However,  as  angle  of  attack  increases,  the  planform  area 
of  one  fin  is  anywhere  from  1  to  4  percent  of  the  body  planform  area  and  it  makes  it 
difficult  to  get  consistently  accurate  values  resulting  from  control  deflection.  This  is 
due  to  the  substraction  process  involved  in  obtaining  the  ACw  due  to  canard  deflec¬ 
tion  and  the  fact  that  there  is  a  certain  amount  of  experimental  error  in  each  of  the 
measurements. 

Second  of  all,  based  on  comparisons  with  data  of  actual  missile  configurations, 
the  effect  of  control  deflection  can  increase  or  decrease  with  angle  of  attack, 
depending  on  the  Mach  number.  At  low  supersonic  and  subsonic  Mach  numbers,  the 
effectiveness  tends  to  get  smaller  with  increasing  angle  of  attack,  whereas  the 
compressibility  effects  at  high  Mach  numbers  cause  the  opposite  effect  to  occur. 
Therefore,  the  smaller  the  effect  from  control  deflection,  the  harder  it  is  to  separate  it 
out  from  that  caused  by  angle  of  attack. 

After  several  weeks  of  reducing  data  and  attempting  to  develop  an  empirical 
model,  it  was  concluded  that  the  inconsistency  of  results  were  such  that  a  different 
approach  for  developing  a  nonlinear  model  of  interference  effects  due  to  control 
deflection  was  desired.  The  approach  taken  was  to  use  {he  AP93  with  the  non- 
linearities  of  wing-alone,  wing-body,  and  body-wing  interference  effects  due  to  angle 
of  attack  included,  use  the  slender-body  estimate;-  of  kw  b  and  kgiw.  for  control 
deflection,  and  derive  empirical  modifications  to  kw  b  based  on  numerical  experi¬ 
ments  compared  to  actual  missile  data.  Because  kw,b  appears  in  the  vortex  lift  on  the 
tail  due  to  canard  or  wing  shed  vortices,  the  numerical  experiments  were  conducted 
with  canard  body-tail  configurations 

Referring  to  Equation  31,  the  vortex  normal-force  coefficient  on  the  tail  is41 

(^N  )  (^N  )  *SviB)  Sin  “  +  F  *Sv(B>  sin8W  i(Sr_rT,AW 
CN  =  - — - — - - -  143) 

™  2„(AR,T(fw-rw)Arf 

Equation  43  has  a  factor  Fthat  multiplies  the  term  due  to  control  deflection  in 
the  wing-tail  vortex  lift.  This  factor  is  needed  in  addition  to  the  nonlinearity  for 
kw(B).  partly  because  the  negative  afterbody  lift  due  to  control  deflection  is  not 
presently  modeled  in  either  AP81  or  AP93.  This  term  is  defined  as 

_  ~4r 
nb<V)  AwV^ 
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where 


_T _ Kw(B)sm«  -f  kW)B)sin8w  /  \ 

V  4  (f  _  r  )  \CN  LAW 

vu0  rw'  l/w 


The  main  reason  this  term  was  not  included  in  the  AP81  code  was  that  it 
required  an  estimate  of  fp,  which  is  the  position  of  the  canard  shed  vortex  at  the  tail. 
Also,  Reference  40  indicated  this  term  was  generally  much  smaller  than  that 
computed  by  Equation  43.  To  account  for  this  term,  a  vortex  tracking  algorithm  is 
needed  or  an  empirical  correction  to  the  term  in  Equation  43.  For  angles  of  attack 
much  greater  than  25  or  30  deg,  a  vortex  tracking  algorithm  may  be  needed. 
However,  up  to  ®  of  about  30  deg,  an  attempt  will  be  made  to  develop  a  nonlinear 
model  of  interference  effects  resulting  from  control  deflection  by  defining  kwiB>  as  a 
function  of  angle  of  attack  and  Mach  number  and  F  as  a  function  of  Mach  number. 


Using  References  42  through  44  and  Reference  40  for  low  Mach  number,  a 
semiempirical  nonlinear  model  for  kw(B>  and  the  parameter  F  were  derived  from 
numerical  experiments.  The  mathematical  model  for  kw(B>  is  based  on  slender-body 
theory  similar  to  kw<B)  and  kB(w>  and  modified  for  angle  of  attack  or  control  deflection. 
In  general,  it  was  found  that 

kW(B)  =  Cl(M)  *Sv<B)  gB  +  C2(i<XW,,M*) 


p  =  c3(m.'v) 


(45) 


More  specifically,  kWiB).  Ci,  C2,  and  F  are  defined  in  Figure  18  for  Mach  numbers 
where  data  are  available.  For  Mach  numbers  less  than  0.8  and  greater  than  4.6,  the 
equations  derived  for  those  conditions  have  been  used.  The  current  method  for  using 
the  empirical  estimate  for  kw<B>  from  Figure  1  is  to  linearly  interpolate  between  Mach 
numbers  for  a  given  value  of  «,  8,  and  Moc. 

In  examining  the  model  in  Figure  18,  it  has  a  lot  of  similarities  to  the  nonlinear 
KwrB)  model  already  discussed;  i.e.,  at  low  angle  of  attack,  slender-body  theory  gives  a 
reasonable  estimate  of  kw(Bi-  However,  as  angle  of  attack  increases,  kw<B>  decreases 
for  low  to  supersonic  Mach  numbers.  For  higher  supersonic  Mach  numbers,  kwiB> 
actually  increases  at  higher  angles  of  attack,  presumably  due  to  compressibility 
effects.  Also,  for  low  angles  of  attack,  a  value  of  Fnear  one  is  found  for  the  vortex  lift 
model,  indicating  again  reasonable  accuracy  using  the  theory  in  Reference  41.  How¬ 
ever,  as  angle  of  attack  is  increased,  F  increases  above  one  for  many  Mach  numbers. 
That  is,  Equation  43  gives  values  of  CnT(Vi  too  small  due  to  control  deflection  of  a 
forward  surface.  As  already  mentioned,  this  is  most  probably  due  to  the  neglect  of  the 
effect  on  the  afterbody  (Equation  44),  which  becomes  more  appreciable  percentage 
wise  compared  to  the  Equation  43  results,  as  angle  of  attack  increases. 
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M<.8 


If  l^w!  £  24.0  — *  kw(B)  —  1 

If  l«wl  >  24.0  — »  kw(g)  =  1.4  [.000794  l«w|2  -  .0933  l«J  4-  2.71] 


F=  1.1 


M=1.1 


If  |«w|  <  1 5.0  — *•  kW(8)  =  1.3[kW(g)]SB 

If  l«wl  >  15.0  -*  kw(B)  =  1.3  (.00087  l«w|2  _  .0825  l*w|  + 1.98] 
F  —  1 . 1 

M  =  1 .5 

If  l«wl  <  10.0 -^kW(B)  =  .9[kW(8)]se 

If  l«wl  >  10.0  -*■  kw(B)  =  9[kw(B)]sB  -  .015{l*wl  -  10.0] 

If  l*wl  £  20.0  F  =  .8 

If  |ocw|  >  20.0  — >  F  =  .8  +  ,10[l«wl  -  20.0] 

M  =  2.0 

If  |ocwl  £  10.0-»kw(B)  =  -9(kw(8)]sB 

If  l*wl  >  10.0  kw(B)  =  •9(kw(B)is8  ~  -005[I<3cwI  -  10.0] 

lfl«wl  £  20.0 -F  =  .8 

If  !<*wl  >  20.0  — *■  F  =  .8  +  -17[l*wl  -  20.0] 

M  =  2.3 

If  l«wl  £  20.0  — *  k^B)  “  -9[kw(B)lsB 

If  l»wl  >  20.0 -^ky^Bj  =  9{kw(B)3s8  ~  005[i*wl  -  20.0] 

If  l«wl  £  30.0  ~  .9 

If  l«wl  >  30.0  — » F  =  .9  +  ,15[l«w|  -  30.0] 

M  =  2.87 

lfl«wl  s  20.0-*kw(B)  =  .9[kw(B)]sB 

If  lawl  20.0  — »  kw(Bj  —  .9[kw(B)]sB  .005(l»wl  —  20.0] 

If  |«J  <  30.0  — F  =  .9 

If  l*wl  >  30.0  -*  F  =  .9  +  .17[|ocwf  -  30.0] 

M  =  3.95 

k  w(B )  =  -8[kw(B)]sB 

If  l«wl  <  40.0  -»  F  =  0.9 

If  |ocwl  >  40.0  -*  F  =  0.9  +  .4(l«wl  -  40.0] 

M>4.6 

If  |ocwl  <  20.0  — >  kw(B)  =  0.75(kw(B)JsB 

If  l«wl  >  20.0— >kw(3)  =  0.75(kw(B)]5B  +  .01[l*wl  -  20.0] 

If  |ocwl  £  35.0  -»  F  =  .9 

If  l«wl  >  35.0  — » F  =  .9  +  ,3{l«wl  -  35.0] 

where  s 
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3.0  SUMMARY  OF  METHODS  IN  AP93 


3.1  AP81  METHODS  USED  IN  AP93 

In  picking  the  methods  to  carry  forward  into  the  AP93  code  from  the  AP81 
version,  several  criteria  were  considered:  (1)  is  the  method  needed  for  performing 
aerodynamic  computations,  (2)  have  the  methods  been  used  by  users  or  has  it  been 
basically  taking  up  computer  storage  and  time  with  little  return  on  the  investment, 
(3)  has  the  AP81  method  been  replaced  by  something  better  in  AP93,  and  (4)  is  the 
accuracy  of  the  method  acceptable  for  most  engineering  purposes.  With  these  criteria 
in  mind,  a  summary  of  the  methods  of  AP81  that  will  be  retained  in  AP93  will  be 
given.  These  methods  have  all  been  well  documented  in  References  2  through  8  so 
the  interested  reader  is  referred  to  those  references  for  details  of  the  methods.  A  brief 
discussion  will  be  given  in  this  report. 

Most  of  the  body-alone  methods  of  AP81  will  remain  in  AP93.  The  nose  wave 
drag  is  currently  computed  at  low  supersonic  Mach  numbers  by  combining  MNT  with 
second-order  Van  Dyke.2,46  This  method,  first  introduced  in  Reference  2,  remains  the 
only  engineering  approach  to  get  accurate  pressure  distributions  and  wave  drag 
estimates  on  blunt  bodies  at  low  supersonic  Mach  numbers  from  Mco  =  12  to  around 
Moo  =  2  to  3  (chosen  by  the  user);  the  SOSET47  48  combined  with  the  improved  MNT 
discussed  earlier  is  used  for  wave  drag  up  to  about  Moo =6.  Above  that  Mach  number, 
the  basic  same  approach  is  used  but  real-gas  and  aeroheating  effects  are  also 
considered,  which  slightly  increases  the  run  time  per  Mach  number.91011  At 
transonic  speeds,  the  wave  drag  is  computed  by  the  method  of  Wu  and  Aoyoma49on 
the  boattail  and  by  a  semiempirical  table  look-up  using  data  generated  by  a  Euler 
code50  for  the  nose.  Of  course,  no  pressure  drag  is  assumed  for  subsonic  flow. 

The  skin-friction  drag  for  both  the  body  and  wings  is  computed  by  the  method  of 
Van  Driest.51  However,  options  for  turbulent  and  laminar  flow  are  offered  as 
discussed  previously.  The  empirical  base  drag  approach  is  retained2  but  replaced  by 
the  improved  method  in  Section  2.0.  The  linearized  lift  and  center-of-pressure 
estimates  developed  for  the  body  alone  will  also  be  retained.  They  are  determined  by 
either  SOSET,47  48  Tsien  first-order  crossflow,52  semiempirical,53  Wu  and  Aoyoma,49 
or  empirical,2  depending  on  the  Mach  number  regions.  The  only  modification  being 
made  to  this  methodology  is  the  center-of-pressure  shift  defined  in  Section  2.0.  The 
nonlinear  normal-force  and  center-of-pressure  methods  of  Allen  and  Perkins34  will  be 
retained  but  extended  as  discussed  in  Section  2.0. 
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Concerning  the  wing  and  interference  aerodynamic  computations,  all  the  static 
aerodynamic  methods  of  References  3  and  7  will  be  retained.  For  normal-force  and 
center-of-pressure  computation,  this  includes  lifting  surface  theory  in  subsonic  flow, 
3-D  thin  wing  theory  (3DTWT)  in  supersonic  flow,54  and  an  option  for  shock 
expansion  strip  theory48  or  3DTWT  in  high  supersonic  or  hypersonic  flow.55-56-57  For 
interference  aerodynamics,  this  includes  slender-body  theory  and  linear  theory,41  and 
line  vortex  theory41  for  wing-tail  interference.  There  is  also  a  term  estimated  on  the 
trailing  edge  of  the  fins  if  they  are  blunt,  which  contributes  to  Fin  base  drag.1 

Concerning  dynamic  derivatives,  only  the  methods  discussed  in  References  2 
(for  the  body  alone)  and  5  (for  the  wing  alone)  will  be  retained.  These  methods 
include  empirical  for  the  body  alone58  and  either  lifting  surface,54  empirical,59  thin 
wing  theory,80'66  or  strip  theory48  for  the  wings.  The  main  reason  for  decreasing  the 
emphasis  on  dynamic  derivatives  is  that  they  do  not  appear  to  be  as  important  to  the 
end  users  of  the  code  as  they  were  in  the  1970’s  when  the  code  was  First  started.  This 
is  because  the  code  is  primarily  used  in  trim  or  particle  aerodynamic  models  where 
static  aerodynamics  and  not  dynamic  derivatives  are  needed.  Also,  guidance  laws 
are  more  robust  today  than  20  years  ago,  and  this  decreases  the  need  for  accurate 
dynamic  derivatives. 


3.2  AP93  OVERALL  METHODS  SUMMARY 

The  methods  used  for  computing  forces  and  moments  in  the  AP93  are 
summarized  in  Tables  6,  7,  and  8.  Note  that  the  code  can  now  be  useful  for  computing 
aerothermal  information  as  well  as  forces  and  moments.  This  means  the  uses  of  the 
code  are  now  for 

•  Providing  inputs  to  flight  dynamics  models  that  estimate  range  or 
miss  distance 

•  Assessing  static  stability  of  various  missile  conFiguration 

•  Assessing  various  design  parameters  in  terms  of  optimizing  the  con¬ 
figuration 

«  Assessing  structural  integrity  using  the  loads  portion  of  the  code 

•  Assessing  aerothermal  aspects  of  a  design  using  the  heat  transfer 
coefficients  at  high  Mach  numbers 

As  seen  in  Tables  6,  7,  and  8,  there  are  many  methods  that  go  into  the  overall  makeup 
of  a  component  buildup  code,  such  as  the  APC.  The  past  20  years  have  shown  that 
this  type  of  code  can  be  quite  useful  when  used  in  preliminary  or  conceptual  design 
studies  to  provide  down  selection  on  many  configuration  alternatives  in  a  fairly 
accurate  and  cost-effective  manner. 
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TABLE  7.  AP93  METHODS  FOR  WING-ALONE  AND  INTERFERENCE  AERODYNAMICS 
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4.0  RESULTS  AND  DISCUSSION 


Some  results  of  component  aerodynamics  or  heat  transfer  coefficients  have 
already  been  presented  during  the  summary  of  the  new  methods  section.  This  section 
of  the  report  will  therefore  focus  on  validation  of  the  overall  AP93  code  in  terms  of 
complete  missile  configurations  where  wind  tunnel  data  are  available.  Where  other 
state-of-the-art  (SOTA)  aerodynamic  codes  have  been  used,  these  resv  lts  will  also  be 
shown  for  comparison  of  the  AP93.  Unfortunately,  funding  will  not  permit 
investigating  another  SOTA  code’s  performance  on  the  example  configurations 
unless  the  results  are  already  available  from  calculations  performed  previously. 

The  first  case  considered  is  the  body-alone  configuration  of  Reference  31.  While 
the  base  pressure  information  was  used  in  the  base-drag  methodology  for  M^2.0,  no 
stability  and  control  information  was  used  in  the  improved  body-alone  method.  The 
database  used  was  primarily  that  of  Reference  36.  Hence,  using  the  AP93  code  on  a 
different  set  of  data  will  help  validate  the  code’s  improved  capability  for  bodies  alone. 
Only  one  case  is  chosen  from  the  many  presented  in  Reference  31.  The  case  selected 
is  one  that  has  a  25-percent  blunt  nose,  is  3.58  calibers  long,  and  has  a  10-caliber 
afterbody.  The  configuration  is  shown  in  Figure  19a. 


TANGENT  OGIVE 


FIGURE  19a.  BODY-ALONE  CONFIGURATION  USED  IN 
VALIDATION  PROCESS^ 
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Results  of  the  comparisons  of  the  AP93,  APH1,  and  experiment  for  the 
configuration  of  Figure  19  are  given  in  Figures  19a  through  19g  in  terms  of  norrnal- 
and  axial-force  coefficients  and  center  of  pressure.  Results  are  given  for  Mach 
numbers  of  0.8,  1.2,  1.5,  2.0,  3.0,  and  4.0.  It  should  he  pointed  out  that  for  Mach 
numbers  of  2.0  and  greater,  only  the  forebody  axial-force  coefficient  is  given.  The 
reason  is  that  no  boundary  layer  trip  was  used  in  any  of  the  tests,  and  the  How  is 
laminar  over  the  entire  body  at  these  Mach  numbers.  This  means  the  base  pressure 
computations  could  have  an  error  built  in  because  of  the  assumption  of  a  turbulent 
boundary  layer  at  the  base.  On  the  other  hand,  for  the  transonic  How  experiments, 
the  Reynolds  number  was  high  enough  to  ensure  turbulent  flow  at  the  base  without  a 
trip.  Thus,  total  axial-force  coefficients  are  shown  in  the  comparison. 

In  viewing  the  comparison  of  AP93  and  AP81  to  the  data  in  Reference  31 .  it  is 
shown  that  AP93  indicates  improvement  over  AP81  at  just  about  all  conditions 
considered  and  at  all  Mach  numbers.  Experimental  data  were  only  available  for 
angles  of  attack  of  about  15  deg,  so  the  theoretical  computations  were  only  shown  up 
to  20  deg.  The  only  area  where  the  AP93  code  appears  to  need  some  improvement  is 
in  the  estimation  of  axial  force  with  angle  of  attack.  The  AP93  code  does  estimate 
base  pressure  as  a  function  of  *  and  has  a  simplified  model  for  skin  friction  as  well  at 
angles  of  attack  up  to  15  deg.  For  supersonic  flow,  the  wave  drag  is  also  computed  as 
a  function  of  a.  However,  it  does  not  take  into  account  separation  regions  on  the  body 
leeside,  which  lower  the  skin  friction  or  even  make  it  negative  on  s./me  portions  of  the 
body,  thus  negating  some  of  the  increased  base  drag  with  angle  of  attack  An 
example  of  this  is  Figure  L9b,  where  the  base  pressure  and  skin  friction  increases 
with  angle  of  attack  are  partially  offset  by  the  flow  separation  on  the  body  leeside, 
which  is  not  accounted  for.  Other  than  this  problem,  the  AP93  code  appears  to  give 
improved  results  for  the  body  alone  over  the  AP81  code. 

A  second  case  considered  is  the  canard-body-tail  case  shown  in  Figure  20a, 
which  was  taken  from  Reference  67.  The  configuration  is  somewhat  of  an  extreme 
case  for  the  body-alone  aerodynamics  because  it  is  100-percent  blunt  and  is  about 
22.3  calibers  long.  The  configuration  tested  in  the  wind  tunnel  has  hangars  attached 
to  the  body  for  aircraft  carry  and  launch.  However,  tests  were  conducted  with  and 
without  the  hangars,  and  the  results  showed  that  CN  and  Cm  were  unchanged  but  Ca 
was  increased.  The  AP93  and  AP81  theoretical  computations  are  compared  to  the 
corrected  data  of  Reference  67,  where  the  hangars  have  been  omitted.  Results  are 
given  in  Figures  20b  through  20s  for  Mach  numbers  of  0.8,  1.2,  2.1,  2.86,  3.95,  and 
4.63  and  at  canard  deflections  of  0,  10,  and  20  deg.  Examining  Figures  20b  through 
20s,  it  is  shown  that  AP93  gives  good  agreement  with  experimental  data  at  just  about 
all  conditions.  Significant  improvements  of  the  AP93  over  the  AP81  are  seen  at  the 
lower  Mach  numbers  and  at  the  higher  Mach  number,  higher  angle-of-attack 
conditions. 
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Anglo  of  Attack  (dog) 


b.  M<*,  =  0.8 


FIGURE  19.  NORMAL- AND  AXIAL- FORCE  COEFFICIENTS  AND 
CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  OF  FIGURE  19a 
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c.  Mcd  =  1.2 

FIGURE  19.  NORMAL-  AND  AXIAL-FORCE  COEFFICIENTS  AND 
CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  OF  FIGURE  19a  (CONTINUED) 
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Angl*  ot  Attack  (deg) 


FIGURE  19.  NORMAL-  AND  AXIAL-FORCE  COEFFICIENTS  AND 
CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  OF  FIGURE  19a  (CONTINUED) 
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Angle  of  Attack  (deg) 


e.  Mr  =  2.0 


FIGURE  19.  NORMAL- AND  AXIAL-FORCE  COEFFICIENTS  AND 
CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  OF  FIGURE  I9a  (CONTINUED) 
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Angle  of  Atlack  (deg) 


f.  M,  =  3.0 


FIGURE  19.  NORMAL- AND  AXIAL-FORCE  COEFFICIENTS  AND 
CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  OF  FIGURE  19a  (CONTINUED) 
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FIGURE  19.  NORMAL-  AND  AXIAL-FORCE  COEFFICIENTS  AND 
CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  OF  FIGURE  19a  (CONTINUED) 
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FIGURE  20a.  CANARD-BODY -TAIL  CONFIGURATION 
USED  IN  VALIDATION  PROCESS1’7 


In  analyzing  why  this  improvement  occurs  at  those  conditions,  it  is  noted  that 
the  aspect  ratio  of  the  tail  surfaces  of  the  configuration  of  Figure  20a  is  about  0.87 
and  that  of  the  canards  is  about  1.7.  Examining  Figure  14  of  Reference  12  (or  Tables 
2  and  3),  it  is  seen  that  the  nonlinearity  in  wing-alone  lift  is  small  for  Mach  numbers 
greater  than  about  1.5.  As  normal  Mach  number  increases,  [M^sin(°c  +  8)]  for  Mach 
numbers  greater  than  about  3.5  to  4.0,  nonlinearity  because  of  compressibility 
becomes  important.  The  bottom  line  is  as  long  as  the  aerodynamics  are  fairly  linear, 
the  AP81  gives  good  results  up  to  moderate  angles  of  attack.  However,  when 
nonlinearities  are  preseal,  tiie  APyj  shows  significant  improvement.  This 
improvement  is  the  greatest  on  the  Figure  20a  configuration  at  low  Mach  number 
because  the  nonlinear  normal-force  term  on  the  canards  is  negative,  whereas  that  of 
the  tails  is  positive.  The  combination  produces  a  strong  couple  in  terms  of  the 
pitching  moment  as  evidenced  by  Figures  20b,  c,  h,  i,  n,  and  o.  A  good  nonlinear 
capability  such  as  that  present  in  the  AP93  is  absolutely  essential  to  get  accurate 
stability  and  control  information  for  these  cases.  Just  examining  Figure  20b,  the 
center  of  pressure  of  the  AP81  at  «  =  20  deg  differs  from  the  experimental  data  by 
-9.4  percent  of  the  body  length  versus  1.3  percent  for  the  AP93. 
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Angia  o<  Attack  (dag) 


e.  M,  =  2.86, 6  =  0  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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Angl*  ot  Attack  (dag) 


Angla  ol  Attack  (Mg) 


f.  Mr  ~3.95,  6-0  DEG 

FIGURE  20.  NORMAL  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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Angla  of  Attack  (dag) 


g.  M,-  =  4.63, 8  =  0  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 


h.  M.  =  0.8. 8  =  10  DEG 

FIGURF  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  r CONTINUED! 
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j.  =  2.1,  f>~  10  DEG 

FIGURE  20  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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Angi*  of  Attack  (d*9) 


k.  M„  =  2.86,  5  =10  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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1.  Mx  =  3.95,  8  =  10  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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10  15  20  25  30 

Angfa  of  Attack  (dag) 


Angla  of  Attack  (dag) 


Angla  of  Attack  (dag) 


m  Mr-  4.63,  fi  =  10  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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FIGURE  20.  NORMAL  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED* 
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Angla  ol  Attack  (dag) 


p.  M^-2.1, 6  =  20 DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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q.  M*  =  2.86, 5  =  20  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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r.  M^  =  3.95,  5  =  20  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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s.  =  4.63,  6  =  20  DEG 

FIGURE  20.  NORMAL-  AND  AXIAL-FORCE  AND  PITCHING 
MOMENT  COEFFICIENTS  FOR  CONFIGURATION 
OF  FIGURE  20a  (CONTINUED) 
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A  third  case  considered  in  the  validation  of  the  development  of  a  nonlinear 
empirical  model  for  kw<Bi  is  a  configuration  tested  at  NASA/LRC42.43  and  is 
representative  of  the  SPARROW  missile.  The  configuration  tested  in  Reference  43  is 
shown  in  Figure  21a.  The  configuration  tested  in  Reference  42  is  just  like  the  one 
tested  in  Reference  43,  except  it  had  wiring  tunnels  and  wave  guides  present.  These 
appendages  add  to  the  normal  force  and  pitching  moment  but  were  not  accounted  for 
in  the  analytical  computations  that  are  presented  in  Figure  21.  The  Reference  43 
configuration  did  not  have  these  appendages  present  and  was  the  main  set  of  data 
used  for  the  nonlinear  empirical  model  for  kw<B>-  These  results  are  distinguished  in 
Figure  21  by  the  fact  that  the  cases  that  had  wave  guides  present  are  indicated. 

Results  of  the  AP81  and  AP93  compared  to  the  experiment  for  the  case  of 
Figure  21a  are  shown  in  Figure  21b  and  c.  Results  are  presented  in  terms  of  C\  and 
Cm  versus  angle  of  attack  for  various  control  deflections  and  Mach  numbers.  As 
shown  in  Figure  19,  the  nonlinear  models  with  and  without  control  defier* 'on  show 
the  AP93  code  agreeing  much  closer  to  the  data  at  all  Mach  numbers  than  the 
linearized  approaches  of  AP81.  On  the  other  hand,  the  fact  that  the  body-alone 
normal  force  of  AP81  had  the  nonlinearities  included  makes  the  comparisons  to 
experimental  data  better  than  it  would  be  otherwise. 

In  examining  Figure  21b,  it  is  seen  that  both  CN  and  Cm  of  AP93  agree  with  the 
experiment  at  8  =  0  and  8  =  10  deg  for  Ma;  =  1.5  and  2.0,  respectively.  C*  and  CM  of  the 
AP81  are  both  considerable  in  error  as  angle  of  attack  increases  above  5  to  10  deg. 
For  M*>  =  2.35,  both  C\  and  Cm  of  AP93  at  8  =  0  and  20  deg  agree  with  the  data. 
Again,  AP81  is  in  considerable  error  at  *  >10  deg,  although  the  error  is  decreasing 
with  increasing  Mach  number.  At  Mx  =  2.87,  both  AP81  and  AP93  give  good  results 
for  Cn  and  Cm  at  8  =  0.  At  8  =  20  deg,  Cn  is  also  predicted  quite  well  by  both  AP93  and 
AP81;  however,  Cm  errors  are  larger  for  AP81  than  AP93.  For  Mx  =  3.95,  AP81  gives 
acceptable  results  for  Cn  and  Cm  up  to  *  =  15  to  20  deg  and  at  both  8  =  0  or  20  deg.  The 
comparison  with  data  gets  worse  above  »  =  20  deg,  whereas  AP93  comparisons  show- 
good  agreement  at  all  *’s  and  8’s.  The  same  statements  basically  hold  true  for  the 
M<x  =  4.6  comparisons. 

Figure  21c  shows  the  comparisons  of  AP81  and  AP93  to  the  data  of 
Reference  42,  which  is  the  same  configuration  as  that  of  Figure  19  except  that  wave 
guides  and  wiring  tunnels  were  attached  to  the  wind  tunnel  model.  As  already 
mentioned,  no  account  was  taken  for  these  appendages  in  the  analytical 
computations.  Note  that  AP93  agrees  much  closer  to  the  data  than  AP81  for  both 
Mx  =  2.3  and  4.6  at  all  values  of  8.  In  comparing  the  wind  tunnel  data  for  the  cases 
with  and  without  appendages,  it  can  be  seen  that  the  appendages  add  only  a  few 
percent  to  the  aerodynamics. 
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FIGURE  2la.  AIR-TO-AIR  MISSILE  CONFIGURATION  USED 
IN  VALIDATION  PROCESS«.« 
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FIGURE  21  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  21a  FOR  VARIOUS  MACH  NUMBERS 
AND  CONTROL  DEFLECTIONS  CONTINUED) 
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1GURE  21  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  2>a  FOR  VARIOUS  MACH  NUMBERS 
AND  CONTROL  DEFLECTIONS 


Mach=4.6,  Delta =-20  Mach=4.6,  Delta=-12 
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FIGURE  21.  NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
CONFIGURATION  OF  FIGURE  21a  FOR  VARIOUS  MACH  NUMBERS 
AND  CONTROL  DEFLECTIONS  (CONTINUED) 


Mach=4,6,  Delta =0  Mach=4.6,  Delta=12 
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A  fourth  case  used  in  the  validation  and  development  of  the  nonlinear  kw,B> 
model  is  shown  in  Figure  2‘2a.  Note  that  in  Figure  22a,  two  configurations  were 
actually  tested44 — one  that  had  a  full-tail  surface  and  a  second  that  had  a  partial 
cutout  removed.  The  APC  will  not  handle  the  partial-wing  configuration  as  it  stands, 
so  an  engineering  model  of  this  wing  must  be  created.  Experience  has  shown  that  the 
lifting  surface  area,  aspect  ratio,  span,  leading  edge  sweep  angle,  and  centroid  of  the 
presented  area  must  be  held  constant.  The  chord  is  varied  so  as  to  meet  these 
constraints.  Hence,  the  configuration  that  represents  the  partial-wing  results  is  the 
body  canard  of  Figure  22a,  plus  the  AP93  representation  of  the  partial  tail  shown  in 
the  lower  right  of  Figure  22a. 


CANARD 


NOSE  RADIUS  1.42 


NOTE:  AU  DIMENSIONS  IN  INCHES.  FUEL  SCALE 
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FIGURE  2‘2a.  CANARD-CONTROLLED  MISSILE 
CONFIGURATION  WITH  FULL-TAIL,PARTIAL-TAIL, 

AND  AP93  REPRESENTATION  OF  PARTIAL  TAIL 
FOR  USE  IN  VALIDATION  PROCESS*1 

Figures  22b  through  22d  present  the  comparison  of  the  AP93  with  wind  tunnel 
test  data.44  Data  were  only  available  at  Moc  —  0.2;  however,  this  compliments  the 
previous  data  set  for  the  SPARROW  missile  in  the  sense  that  no  subsonic  data  were 
available  for  that  case.  The  full-tail  and  partial-tail  results  are  denoted  on  the  figure. 
Some  results  were  available  from  Reference  44  for  the  Missile  Datcom. *5  These 
results  are  also  shown  where  available.  As  seen  in  the  figure,  the  AP93  gives 
improved  results  for  pitching  moment  and  normal  force  for  most  conditions  compared 
to  the  Missile  Datcom.  While  center  of  pressure  is  not  shown,  the  AP93  computations 
are  generally  within  the  goal  of  ±4  percent  of  the  body  length.  For  example,  at 
'x  =30  deg,  8  = -20  deg,  xcp  for  the  data,  AP93  and  Missile  Datcom  are  5.39,  4.91,  and 
3.75  calibers,  respectively,  with  respect  to  the  moment  reference  point.  This 
represents  errors  of  2.1  and  7.3  percent  of  the  body  length,  respectively,  for  the  AP93 
and  Missile  Datcom  codes. 
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FIGURE  22b.  COMPARISON  OF  AP93  TO  WIND  TUNNEL  DATA  AND 
MISSILE  DATCOM  FOR  NORMAL-FORCE  AND  PITCHING  MOMENT 
COEFFICIENTS  OF  FIGURE  22a  CONFIGURATION 


FULL  TAIL  „  PARTIAL  TAIL 
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FIGURE  22d.  COMPARISON  OF  AP93TO  WIND  TUNNEL  DATA  AND 
MISSILE  DATCOM  FOR  NORMAL  FORCE  AND  PITCHING  MOMENT 
COEFFICIENTS  OF  FIGURE  22a  CONFIGURATION  (CONTINUED) 
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Two  body-tail  cases  were  considered  in  the  validation  process.  The  first  of  these 
is  shown  in  Figure  23a  and  is  taken  from  Reference  68.  Several  wing  planforms  were 
tested  in  Reference  68,  but  only  one  example  will  be  shown.  The  body  was  just  over 
13  calibers  long  with  a  0.54-caliber,  4-deg  boattail  at  the  rear.  For  these  tests,  a 
boundary  layer  trip  was  used,  and  data  were  obtained  at  *  up  to  25  deg  at  Mach 
numbers  1.6,  2.36,  and  2.86  and  for  8  =  0  and  10  deg. 


MOMENT  REFERENCE  CENTER 


FIGURE  ’23a.  BODY-TAIL  CONFIGURATION  USED  IN  VALIDATION  PROCESS'* 

(DIMENSIONS  IN  INCHES) 


Figures  23b  through  23c  present  the  wind  tunnel  measurements,  AP93  and 
AP81  for  the  Mach  1.6  and  2.86  cases  at  control  deflections  of  0  and  10  deg.  Results 
are  presented  in  terms  of  normal,  forebody  axial,  and  pitching  moment  coefficients  as 
a  function  of  angle  of  attack.  The  only  real  improvement  of  the  AP93  code  compared 
to  the  AP81  version  is  in  pitching  moment  at  high  <*  and  at  transonic  speeds.  The 
reason  for  the  slight  improvement  in  the  configuration  is  a  body-tail  case  with  a 
moderate  aspect  ratio  tail  surface.  This  type  of  configuration  is  typical  of 
configurations  where  the  AP81  code  performed  satisfactorily. 
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Angt«  of  Attack  (<3tg) 

c.  M,  =  1.6, 8  =  10  DEG 


FIGURE  23.  AERODYNAMICS  OF  A  BODY-TAJI-  CONFIGURATION 
USED  IN  VALIDATION  PROCESS'*  (CONTINUED) 
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of  Attack  (d«g) 

d.  =  2.86, 8  =  0  DEG 


FIGURE  23.  AERODYNAMICS  OF  A  BODY-TAIL  CONFIGURATION 
USED  IN  VALIDATION  PROCESS-’8  (CONTINUED) 


Angla  ot  Attack  (dag) 


e.  M,  =  2.86,8  =  ioi)EG 

FIGURE  23.  AERODYNAMICS  OF  A  BODY-TAIL  CONFIGURATION 
USED  IN  VALIDATION  PROCESS'^  (CONTINUED) 
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The  second  body-tail  configuration  is  shown  in  Figure  24a,  and  the  wind  tunnel 
test  data  along  with  details  of  the  test  are  given  in  Reference  69.  The  configuration  of 
Figure  24a  is  10  calibers  long  with  an  aspect  ratio  tail  of  about  1.0.  Again,  this 
configuration  is  typical  of  cases  where  the  AP81  code  has  performed  well  in  the  past. 
Results  in  terms  of  CN  and  CM  as  a  function  of  *  at  Me*  =  0.6,  0.9,  and  1.3  are  shown  in 
Figure  24b  for  the  AP93,  AP81,  and  data.  The  AP93  shows  some  improvement  over 
the  AP81  for  pitching  moment  at  high  *,  but  normal  force  is  about  the  same.  The 
reason  for  this  improvement  is  the  nonlinear  lift  loss  methodology  included  in  the 
AP93  that  was  not  available  in  the  AP81  code. 
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FIGURE  24a.  BODY-TAIL  CONFIGURATION  USED  IN  VALIDATION  PROCESS^ 
(DIMENSIONS  IN  INCHES,  FULL  SCALE) 
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The  last  two  configurations  used  for  the  validation  process  are  from 
Reference  70.  These  configurations  feature  somewhat  unconventional  missile  fin 
arrangements  that  should  help  in  validating  the  AP93.  The  first  of  these  con¬ 
figurations  is  shown  in  Figure  25a  and  consists  of  a  fairly  large  wing  followed  by  a 
much  smaller  tail  that  is  used  for  control.  The  ratio  of  the  wing  and  tail  planform 
areas  to  that  of  the  body  is  16.25  and  2.99,  respectively.  The  aspect  ratios  of  these 
lifting  surfaces  are  0.92  and  3.63,  respectively,  and  their  r/s  values  are  0.23  and  0.26. 
Example  computations  are  shown  for  Mach  numbers  1.61,  2.0,  2.86,  3.95,  and  4.63 
and  at  tail  deflections  of  0,  -10  deg,  and  -20  deg. 

Comparisons  are  made  with  AP93,  AP81,  and  the  Reference  70  wind  tunnel 
data  for  lift,  drag,  and  pitching  moment  coefficients  in  Figures  25b  through  25p.  For 
this  missile  configuration,  the  AP93  shows  slight  improvement  in  predictions  over 
AP81.  Both  codes  give  fair  to  good  engineering  estimates  of  aerodynamics  at  all 
conditions  considered.  The  area  where  the  AP93  code  still  needs  slight  improvement 
is  in  the  Mach  2.0  to  3.0  range  at  higher  angles  of  attack  and  for  AR  of  about  one. 

In  attempting  to  analyze  why  this  problem  occurs,  the  wing-alone  data  bases  of 
References  39  and  40  were  reviewed  again.  The  nonlinear  wing-lift  model  of  AP93  is 
based  on  the  Reference  39  data  for  Mach  numbers  of  1.6  and  greater  and  on 
Reference  40  data  for  Mach  numbers  of  1.2  and  lower.  Reference  39  data  were  com¬ 
pared  with  that  of  other  data  sets,  including  Reference  40.  At  Mach  2  and  AR=  1.0, 
the  Reference  39  data  are  up  to  15  percent  lower  than  the  Reference  40  wing-alone 
data.  It  is  suspected  that  since  the  Reference  39  forces  and  moments  were  not 
measured  but  integrated  from  pressure  data,  and  since  the  wings  were  required  to  be 
fairly  thick  to  house  the  pressure  taps,  that  this  is  the  reason  for  the  difference.  It  is 
also  suspected  that  around  M  =  1.5  to  2.5,  the  loss  in  lift  due  to  thickness  is  larger 
than  at  other  Mach  numbers.  In  hindsight,  it  would  have  been  better  to  develop  the 
wing-alone  nonlinear  model  based  on  an  average  of  the  two  data  bases.  In  future 
work,  this  will  be  done. 
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FIGURE  25a.  WING  BODY  TAIL  CONFIGURATION'  USED 
IN  VALIDATION  PROCESS™ (DIMENSIONS  IN  INCHES) 
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FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS™ 


4-52 


Angle  ol  Attack  (d«g) 


c.  M^  2.0, 5  =  0  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70  (CONTINUED) 
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d.  M  ,  =  2.86, 8  =  0  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS7*)  (CONTINUED) 
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e.  M,  =  3.95,8  =  0  DEG 


FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING  BODY  TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70 (CONTINUED) 
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f.  Mr- 4.63, 5  =  0  DEG 


FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY  TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70 (CONTINUED) 
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FIGURE  25,  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS?*) (CONTINUED) 


4-57 


FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING  BODY  TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS^ (CONTINUED) 


4-58 


ot  Attack  (d#^) 


i.  M„  =  2.86, 6-10  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70 (CONTINUED) 


Angl*  ol  Attach  (4*g) 


j.  -3.95,  8  =  -10  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70  (CONTINUED) 


An$«  of  Attack  (d«cj) 


k.  M„  =  4.63,  6=10  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70 (CONTINUED) 


Angle  ot  Attack  (deg) 


1.  M„  =  1.61, 8= -20  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS™  (CONTINUED) 


Anglt  o'  Alack 


m.  =  2.0,  fi  =  -20  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS7*)  (CONTINUED) 


Angia  ot  Attack  (dag) 


n.  M,  =  2.86,5  =  -20  DEG 

FIGURE  25.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70  (CONTINUED) 
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The  final  configuration  considered  in  the  validation  process  is  the  canard-wing- 
body  configuration  shown  in  Figure  26a.  As  mentioned  earlier,  the  wind  tunnel 
results  for  this  configuration  were  also  given  in  Reference  70.  Again,  results  are 
shown  for  Mae  =  1.61,  2.0,  2.86,  3.95,  and  4.63  at  control  deflections  of  0-,  10-,  and 
20-deg.  Comparisons  of  the  AP93  to  the  AP81  and  experiment  show  the  AP93  to  he 
slightly  superior  to  the  AP81  at  most  conditions  considered  for  the  forces  and 
moments.  The  greatest  improvement  is  shown  for  the  low  Mach  number  (1.61) 
pitching  moment  at  all  control  deflections.  The  AP81  actually  is  better  than  the 
AP93  at  M,-2.0,  where  the  linear  theory  actually  out  performs  the  nonlinear 
theory.  Again,  it  is  suspected  that  the  previous  discussion  for  the  Figure  25a  con¬ 
figuration  also  applies  here.  Since  this  has  been  fairly  consistently  true  for 
configurations  that  have  large  lifting  surfaces  with  low  values  of  r/s,  some  adjust¬ 
ment  in  the  nonlinear  wing-alone  lift  methodology  is  still  needed. 

Reference  12  considered  several  cases  that  had  low  aspect  ratio  lifting  surfaces. 
For  configurations  such  as  those,  the  AP93  shows  substantial  improvement  in  the 
aerodynamic  estimation  process.  Most  of  the  data  for  these  type  of  configurations  is 
either  classified  or  limited  distribution;  therefore,  no  additional  results  are  shown. 
However,  the  conclusions  in  that  report  still  hold  true. 


FIGURE  2fia.  CANARD  WING  BODY  CONFIGURATION  USED  IN  VALIDATION  PROCESS'" 

(DIMENSIONS  IN  INCHES) 


4-67 


NSWCDD/TR-93/91 


b.  1.61, 6-0  DEG 


FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70 
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FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS7*1  (CONTINUED) 
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f.  M.,  =4.63, 8  =  0  DEG 


FIGURE  26.  LIFT.  DRAG,  AND  PITCHING  MOMENT  COEFFIC  IENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS7'* ( CONTINUED) 
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Angf«  o I  Attack  (<J«g) 


g.  M«=  1.61,6  =  +10  DEG 


FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING  BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70  (CONTINUED) 
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Angl«  of  Artadi  (d»g) 


h.  M*  =  2.0,  8=  +  10  DEG 


FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS7®  (CONTINUED) 
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FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70  (CONTINUED) 
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FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD- WING- BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS-0  (CONTINUED) 
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Angl«  ol  Attack  (a«g) 


k.  M„=4.63,  8=  10  DEG 

FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70 (CONTINUED) 
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Angi#  of  Attack  (dag) 


1.  M*  =  1.61, 8=20  DEG 


FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS^ (CONTINUED) 
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Angle  of  Attack  (deg) 

m.  Mj.  =  2.0,  6  =  20  DEG 


FIGURE  26.  LIFT,  DRAG,  AND  PITCHING  MOMENT  COEFFICIENTS 
OF  A  CANARD- WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS™  (CONTINUED) 
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OF  A  CANARD-WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS70  (CONTINUED) 
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5.0  SUMMARY  AND  RECOMMENDATIONS 


In  summary,  an  improved  APC  (AP93)  has  been  developed  This  report 
summarizes  the  methodology  for  the  new  technology  previously  developed  and 
documented,  presents  in  more  detail  an  undocumented  nonlinear  method  for  wing- 
body  interference  due  to  control  deflection,  and  compares  the  new  AP93  code  to  the 
former  APC  (AP81)  and  experimental  data. 

New  technology  developed  and  summarized  includes  extension  of  the  SOSET  to 
include  real-gas  effects;  a  new  engineering  method  to  calculate  heat  transfer  coeffi¬ 
cients  and  adiabatic  wall  temperature;  new  approximate  pressure  prediction  schemes 
for  blunt  and  sharp  bodies;  an  improved  body-alone  nonlinear  lift  prediction  method; 
new  nonlinear  lift  prediction  methods  for  wing-alone,  wing-body  interference,  body¬ 
wing  interference,  and  wing-body  interference  due  to  control  deflection;  and  a  new 
base  drag  database  to  estimate  the  effects  of  angle  of  attack,  fin  location,  and  size, 
along  with  an  improved  empirical  prediction  model. 

New  capabilities  of  the  code  include  the  ability  to  calculate  information  to  be 
used  for  conducting  engineering  estimates  of  heat  transfer  at  high  Mach  numbers, 
the  ability  to  use  the  code  to  get  accurate  estimates  of  nonlinear  aerodynamics  for  low 
aspect  ratio  and  low  or  high  Mach  number  aerodynamics  up  to  angles  of  attack  of 
30  deg,  and  improved  axial-force  estimations. 

Comparison  of  the  AP93  code  to  the  AP8I  version  and  experimental  data 
showed  the  AP93  code  on  average  reduces  the  errors  in  normal-force  and  center-of- 
pressure  estimation  considerably  for  most  configurations  and  gives  slightly  improved 
estimates  of  axial-force  coefficient.  Comparing  the  AP93  code  to  available  computa¬ 
tions  from  the  external  literature  of  other  SOTA  codes  showed  the  AP93  to  be  as  good 
as  or  superior  to  other  codes  for  planar  aerodynamics.  Regions  of  superiority  of  the 
AP93  code,  in  general,  were  the  same  weak  areas  of  the  AP81  code;  i.e.,  lower  aspect 
ratio,  lower  Mach  number,  higher  angle  of  attack,  pitching  moments,  and  hypersonic 
Mach  numbers. 

The  computational  time  of  the  AP93  is  about  the  same  as  the  AP81  version 
(about  1  sec  per  case  on  a  CDC  995  computer).  This  is  because  many  of  the  unused  or 
seldom  used  techniques  of  the  AP81  version  were  eliminated  and  replaced  with  more 
recent  technology.  There  are  a  few  new  inputs,  but  these  are  minimal.  In  general, 
the  AP93  has  maintained  the  same  level  of  computational  time  and  ease  of  use  (less 
than  one-half  day  for  an  experienced  user  to  set  up  a  configuration  for  computational 
purposes)  but  is  much  more  accurate  and  robust  than  the  AP81  version.  The  code  will 
be  made  available  to  legitimate  requesting  users  at  no  charge  to  them. 
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While  the  AP93  shows  substantial  improvement  over  the  AP81  code  for  most 
static  aerodynamics  on  most  configurations,  there  are  still  areas  needing  improve¬ 
ment.  These  include 

•  Refinements  in  the  wing-alone  lif4  curves  for  wings  that  have 
moderate  aspect  ratios 

•  Refinements  in  the  nonlinear  wing-body  interference  lift  at 
moderate  supersonic  Mach  numbers 

•  Extension  of  the  methodology  to  higher  angle  of  attack 

•  Extension  of  the  methodology  to  the  0  =45-deg  plane  in  addition  to 
the  0  =0-deg  plane 

These  improvements  will  be  added  as  funding  permits. 
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7.0  SYMBOLS  AND  DEFINITIONS 


Ap 

Planform  area  of  the  body  or  wing  in  the  crossflow  plane  (ft2) 

Aref 

Reference  area  (maximum  cross-sectional  area  of  body  if  a  body  is 
present  or  planform  area  of  wing  if  wing  alone)  (ft2) 

Aw 

Planform  area  of  wing  in  crossflow  plane  (ft2) 

a 

Speed  of  sound  (ft/sec) 

AR 

Aspect  ratio  =  b2/Aw 

b 

Wing  span  (not  including  body)  (ft) 

CA 

Axial  force  coefficient 

Cm 

Pitching  moment  coefficient  (based  on  reference  area  and  body 
diameter  if  body  present  or  mean  aerodynamic  chord  if  wing 
alone) 

Cn 

(  Normal  Force  \ 

Normal  force  coefficient  1  1 

v  1/2  p  V2A  ./ 

rx  ao  ref 

Cn^vi 

Negative  afterbody  normal-force  coefficient  due  to  canard  or  wing 
shed  vortices 

Additional  normal-force  coefficient  on  body  due  to  presence  of 
wing 

ACNb,w, 

Additional  normal-force  coefficient  on  body  due  to  a  control 
deflection  of  the  wing 

CNl 

Linear  component  of  normal-force  coefficient 

Cnn, 

Nonlinear  component  of  normal-force  coefficient 

Negative  normal-force  coefficient  component  on  tail  due  to  wing 
or  canard  shed  vortex 
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CN 

ACNw(U, 

Cn, 

Cp 


CpB 

^Pb^,S,Vc,X/C 

CP0 


Cr 


Ct 

d 

dref 

e 

F 


Fi,  Fa,  F3 


fw»  fit 


H 

Ho 


Normal-force  coefficient  of  wing  in  presence  of  body 

Additional  normal-force  coefficient  of  wing  in  presence  of  body 
due  to  a  wing  deflection 

Normal-force  coefficient  derivative 

Pressure  coefficient  (  — — \ 

'  1/2  p  V2  ' 

Base  pressure  coefficient 

Base  pressure  coefficient  with  no  fins  present  and  at  angle  of 
attack 

Base  pressure  coefficient  with  fins  present  of  some  t/c,  x/c,  6,  and 
body  at  some  « 

Stagnation  pressure  coefficient 
Root  chord  (ft) 

Tip  chord  (ft' 

Body  diameter  (ft) 

Reference  body  diameter  (ft) 

Internal  energy  (ffc2/sec2) 

Dimensionless  empirical  factor  used  in  tail  normal-force 
coefficient  term  due  to  canard  or  wing  shed  vortices  to 
approximate  nonlinear  effects  due  to  a  control  deflection 

Symbols  defining  parameters  used  in  base  drag  empirical  model 

Lateral  location  of  wing  or  tail  vortex  (measured  in  feet  from  body 
center  line) 

Heat  transfer  coefficient  based  on  wall  local  temperature 
(ft-lb)/(ft2-sec-QR) 

Total  enthalpy  (ft2/sec2) 
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Heat  transfer  coefficient  based  on  wall  local  specific  enthalpy 
[slug/(ft2-sec)l 

Specific  enthalpy  (ft2/sec2) 

Adiabatic  wall  specific  enthalpy  (ft2/sec2) 

Specific  enthalpy  at  outer  edge  of  boundary  layer  ( ft2/sec2) 

Height  of  wing  or  canard  shed  vortex  at  tail  center  of  pressure  (ft) 

Specific  enthalpy  at  wall  (ft2/sec2) 

Reference  value  of  specific  enthalpy  (ft2/sec2)  (see  Equation  16) 
Tail  interference  factor 

Empirical  factor  defined  in  wing-alone  nonlinear  normal-force 
coefficient  term 

Ratio  of  additional  body  normal-force  coefficient  derivative  due  to 
presence  of  wing  to  wing-alone  normal-force  coefficient  derivative 
at  8  =  0  deg 

Ratio  of  normal-force  coefficient  derivative  of  wing  in  presence  of 
body  to  that  of  wing  alone  at  8  =  0  deg 

Ratio  of  additional  body  normal-force  coefficient  derivative  due  to 
presence  of  wing  at  a  control  deflection  to  that  of  the  wing  alone 
at  *  =  0 


Ratio  of  wing  normal-force  coefficient  derivati  i  in  presence  of 
body  due  to  a  control  deflection  to  that  of  wing  alone  at «  ^  0  deg 

Value  of  kw(B)  calculated  by  slender-body  theory  at «  =0 

Nonlinear  corrections  to  Kb(w>  and  Kwib>  due  to  higher  angle  of 
attack 

Length  (ft) 

Nose  length  (can  be  in  calibers  or  feet) 

Linear  theory 
Mach  number  =  V/a 
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Mn 

Ni,Nt 


P 

Pc 

Pr 

qw 


Qw,],  9w,t 


R 


Re 

(Re)c 


ReD 


r 


rn 


rw>  ft 


r/s 

s 


SB 

T 

Taw,  To,  Tw 


t/cr 


Normal  Mach  number  to  body  axis  =  M  sin  « 

Transformation  factors  used  in  Eckert  reference  enthalpy  to 
approximate  three-dimensional  effects  for  laminar  and  turbulent 
flow  ( =  3  and  2,  respectively) 

Pressure  (lb/ft.2) 

Pressure  of  a  cone  of  given  half  angle  (lb/ff2) 

Prandtl  number 

Heat  transfer  rate  (ft-ib)/(ft2-sec)  at  wall 

Heat  transfer  rate  at  wall  for  laminar  or  turbulent  flow, 
respectively 

Gas  constant  [for  air  R=  1716  ft-lb/(slug-°R)| 

PV1 

Reynolds  number  = 

Critical  Reynolds  number  where  flow  transitions  from  laminar  to 
turbulent  flow 

Reynolds  number  based  on  diameter  of  wing  leading  edge 
bluntness 

Radius  of  body  (ft) 

Radius  of  nose  tip  (ft) 

Radius  of  body  at  wing  or  tail  locations 

Ratio  of  body  radius  to  wing  or  tail  semispan  plus  the  body  radius 

Distance  along  body  surface  in  SOSET  (also  wing  or  tail  semispan 
plus  the  body  radius  in  wing-body  lift  methodology) 

Slender-body  theory 

Temperature  (°R) 

Adiabatic  wall,  total,  and  wall  temperature,  respectively 
Tail  thickness  to  its  root  chord 
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t/d 

V 

Ve 

Vp 

x 

x/c 


xcp 

XL,  XT 

Z 

a 

ac 

aD 

aW>  aT 

5 

8eq 
8W,  8t 


Tail  thickness  to  body  diameter 
Velocity  (ft/sec) 

Velocity  at  edge  of  boundary  layer  (ft/sec) 

Velocity  parallel  to  leading  edge  of  wing  (ft/sec) 

Distance  along  the  axis  of  symmetry  measured  positive  aft  of  nose 
tip  (feet  or  calibers) 

Parameter  used  in  base  drag  methodology  to  represent  the 
number  of  chord  lengths  from  the  base  (measured  positive 
upstream  of  base) 

Center  of  pressure  (in  feet  or  calibers  from  some  reference  point 
that  can  be  specified) 

Laminar  and  turbulent  flow  lengths  on  body  (ft) 

Compressibility  factor 
Angle  of  attack  (degrees) 

Angle  of  attack  where  wing-body  interference  factor  starts 
decreasing  from  its  slender-body  theory  value  (degrees) 

Angle  of  attack  where  the  wing-body  interference  factor  reaches  a 
minimum  (degrees) 

Local  angle  of  attack  of  wing  or  tail  («  +  8W  or  «  +  5T,  respectively, 
in  degrees) 

Control  deflection  degrees 

Angle  between  a  tangent  to  the  body  surface  at  a  given  point  and 
the  velocity  vector  (degrees) 

Deflection  of  wing  or  tail  surfaces  (degrees),  positive  leading  edge 
up 


4>  Circumferential  position  around  body  where  $  =  0  is  leeward 

plane  (degrees) 

A  Taper  ratio  of  a  lifting  surface  =  ct/cr 
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n 


no 

Ho,  p* 


P,  Po,  p* 


Y 

0 

A 

00 

2- D 

3- D 

3DTWT 

AP81 

AP93 

APC 

GSET 

IMNT 

MNT 

NASA7LRC 

NSWCDD 

SE 

SOSET 

SOTA 


Parameter  used  in  SOSET  and  also  used  in  viscous  crossflow 
theory  for  nonlinear  body  normal  force  (in  this  context,  it  is  the 
normal  force  of  a  circular  cylinder  of  given  length-to-diameter 
ratio  to  that  of  a  cylinder  of  infinite  length) 

Value.of  q  in  viscous  crossflow  theory  for  =  0 

Viscosity  coefficient  at  stagnation  or  reference  conditions, 
respectively  (slug/ft-sec) 

Density  of  air  at  local,  stagnation,  or  reference  conditions, 
respectively  (slugs/ft-'fi 

Specific  heat  ratio 

Local  body  slope  at  a  given  point  (degrees) 

Leading  edge  sweep  angle  of  wing  or  tail  (degrees) 

Free-stream  conditions 
Two  dimensional 
Three  dimensional 
3-D  thin  wing  theory 
Aeroprediction  1981 
Aeroprediction  1993 
Aeroprediction  code 
Generalized  shock-expansion  theory 
Improved  modified  Newtonian  theory 
Modified  Newtonian  theory 

National  Aeronautics  and  Space  Administration/Langley  Re¬ 
search  Center 

Naval  Surface  Warfare  Center,  Dahlgren  Division 
Shock  expansion 

Second-order  shock-expansion  theory 
State  of  the  art 
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